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ABSTRACT 

We present the discovery of KELT- lb, the first transiting low-mass companion from the wide-field Kilo- 
degree Extremely Little Telescope-North (KELT-North) transit survey, which surveys ~ 40% of the northern 
sky to search for transiting planets around bright stars. The initial transit signal was robustly identified in the 
KELT-North survey data, and the low-mass nature of the occultor was confirmed via a combination of follow- 
up photometry, high-resolution spectroscopy, and radial velocity measurements. False positives are disfavored 
by the achromaticity of the primary transits in several bands, a lack of evidence for a secondary eclipse, and 
insignificant bisector variations. A joint analysis of the spectroscopic, radial velocity, and photometric data 
indicates that the V = 10.7 primary is a mildly evolved mid-F star with r e g = 6518 ±50 K, logg* = 4.229^Q glg 
and [Fe/H] = 0.008 ± 0.073, with an inferred mass M* = 1 .324 ± 0.026 M and radius = 1 .462!^ R Q . 
The companion is a low-mass brown dwarf or a super-massive planet with mass Mp = 27.23^ 4^ Mj up , radius 
Rp = 1.1 10^022 ^Jup, surface gravity loggp = 4.738^0 Q23> and a density pp = 24.7*1'% g cm" 3 . The companion is 
on a very short (~ 29 hour) period circular orbit, with an ephemeris ^.(BJDtdb) = 2455909.292797 ± 0.00024 
and P = 1 .2175007 ±0.000018 d, and a semimajor axis of a = 0.02466 ±0.00016AU. KELT- lb receives a large 
amount of stellar insolation, with (F) = 7.81^33 x 10 9 erg s" 1 cm" 2 , implying an equilibrium temperature as- 
suming zero albedo and perfect redistribution of T eq = 2422+?,;? K. Upper limits on the secondary eclipse depth 
in i and z bands indicate that either the companion must have a non-zero albedo, or it must experience some 
energy redistribution. Comparison with standard evolutionary models for brown dwarfs suggests that the radius 
of KELT- lb is likely to be significantly inflated. Adaptive optics imaging reveals a candidate stellar companion 
to KELT-1 with a separation of 588 ± lmas, which is consistent with an M dwarf if it is at the same distance as 
the primary. Rossiter-McLaughlin measurements during transit imply a projected spin-orbit alignment angle 
A = 2 ± 16 degrees, consistent with the orbit pole of KELT-lb being aligned with the spin axis of the primary. 
Finally, the vsin/» = 55.4 ±2.0 km s" 1 of the primary is consistent at ~ 2 a with tidal synchronization. Given 
the extreme parameters of the KELT- 1 system, we expect it to provide an important testbed for theories of the 
emplacement and evolution of short-period companions, as well as theories of tidal dissipation and irradiated 
brown dwarf atmospheres. 

Subject headings: planetary systems, stars: individual: KELT-1, TYC 2785-2130-1, techniques: photometric, 
techniques: spectroscopic 
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1. introduction 

The most information-rich exoplanetary systems are those 
in which the companion happens to transit in front of its par- 
ent star. Transiting systems are enormously useful for en- 
abling detailed measurements of a seemingly endless array of 
physical propert ies of extrasolar pl anets and their host stars 
(see reviews by IWinnl 120091 12010I) . The most basic prop- 
erties that can be measured using transiting planets are the 
planet mass and radius, and so average density. These pa- 
rameters alone allow for interesting constrai nts on the in- 
ternal composition and structure of planets (jGuillotJ 120051 
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20111) . In addition to these basic parameters, transiting plan- 
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ets en a ble the study of the i r atmo s pheres (|Seager & Sasselovl 
~~ aD l200l IVidal-M adiaret'all " 
>n ([Demin 

20051: ICharbonneau et all 120051: |KnutsonetaD[2008). 
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and thermal emission (Demin g et al.l 

They 



also allow measurement of planetary and stellar o blate- 
ness, r otation rate, and sp i n-orbit alignment dSeager & Hui 
' 20021: ISpiegel et all 12007b ICarter & Wind 1201 Ot lRossitei 
19241 iMcLaughlinl 1 19241: Iwinn et al.l 12005b Gaudi & Wine 
20061) . Transiting pla nets may also be searched for associ 



ated rings and moons dBrown et al.l 120011: iBarnes & Fortnevl 
120041: iTusnski & Valid 1201 lb . Further, variations in transit 
timing may indicate the pr e sence of other bodies in the system 
(Holman & Murravl 120051: lAgol et al.l l2005t ISteffen & Agoll 
2005 [ iFord & Gaudil 120061: iFord & Holmanl 12007b iKipping 
2009). With sufficiently precise observations, one may con- 
strain the presence of planets wi t h masses smaller than that of 
the Earth (lAgol & Steffenll2006l: ICarter & WinnllMToh . 

The high scientific value of transiting planet systems mo- 
tivated the first dedicated wide-field transit surveys, which 
by now have ident ified over 100 transiting syste ms (TrES, 
Alonso et al J l2004t XO, iMcCullough et all 120061: HATNet 
Bakos et al.l2007b SuperW ASP jCollier Cameron et al.l2007al 
PES. lAlsubai et al.ll20llh . Although there is substantial di- 
versity in their design, strategy, and sensitivity, these surveys 
can be grossly characterized as having relatively small cam- 
eras with apertures of order 10 cm, and detectors with rela- 
tively wide fields-of-view of tens of square degrees. These 
surveys are primarily sensitive to giant, close-in planets with 
radii Rp > 0.5/?j up and periods of P < 10 d, orbiting relatively 
bright FGK stars with V ~ 1 - 1 2 . 

Th e space-based missi ons C0R0T (Baglin]|200i and Ke- 
pler dBorucki et al.l l2010) have dramatically expanded the pa- 
rameter space of transit surveys, enabling the detection of 
planets with sizes down to that of the Earth and below, plan- 
ets with periods of several years, and planets orbiting a much 
broader range of host stars. Furthermore, their large tar- 
get samples have allowed the detection of rare and there- 
fore interesting planetary systems. These missions have al- 
ready announced over 50 confirmed planets, and the Ke- 
pler mission has ann ounced an additional ~ 2300 candidates 
(Batal ha et alJ 12012) . most of which are smaller than Nep- 
tune. Notable individual discoveri es include the fi rst detec- 
tion of a transiting Super-Earth (Leger et al1 l2009l) . the de- 
tection of a 'temper ate' gas giant wit h a relatively long pe- 
riod of — 100 days dDeeg et all 1201 Ol) the first multi-planet 
transiting systems ( Steff en et all 1201 Ot iHolman et all l2010t 
lLatham et al] 1201 lb iL issauer et al. 201 1 J), the first circumbi- 
nary planets (iDoyle et al.l 1201 lb IWelsh et"aT1 20121) . and the 



detection of planets wit h radius of < R§ (Muirhea d et al.l 
|20T1 iFressin eUd1l20l2l) . 

Although Kepler and C0R0T have revolutionized our un- 
derstanding of the demographics of planets, the opportunities 
for follow-up of the systems detected by these missions are 
limited. By design, both missions primarily monitor relatively 
faint stars with V > 12. Consequently, many of the follow-up 
observations discussed above that are generically enabled by 
transiting systems are not feasible for the systems detected 
by Kepler and C0R0T. Detailed characterization of the major- 
ity of these systems will therefore be difficult or impossible. 
There is thus an ongoing need to discover transiting planets 
orbiting the bright stars, as well as to increase the diversity of 
such systems. 

All else being equal, the brightest stars hosting transiting 
planets are the most valuable. Larger photon flux permits 



more instruments and/or facilities to be employed for follow- 
up, allows subtler effects to be probed, reduces statistical un- 
certainties, and generally allows for improved or more exten- 
sive calibration procedures that help to control systematic er- 
rors. Furthermore, brighter stars are also easier to charac ter- 
ize, and are more likely to have pre-existing information, such 
as proper motions, parallaxes, metallicities, effective temper- 
atures, angular diameters, and broadband colors. 

The majority of the brightest (V < 8) FGK dwarfs in the 
sky have been monitored using precision radial velocity sur- 
veys for many years, and as a result most of the giant planets 
with periods of less than a fe w years orbiting thes e stars have 
already been discovered (e.g., Wri ght et al.l l2012). A smaller 
subset of these stars have been monitored over a shorter time 
baseline with the sensitivity needed to detect Neptune- and 
SuperEarth-mass planets. Because of the low a priori transit 
probability for all but short period planets, the transiting sys- 
tems constitute a very small fraction of this sample. To date, 
seven planets first discovered via radial velocity have subse- 
quently been discovered to also transit; all of the host stars 
for these planets are brighter than V = 9. Although there ar e 
projects that aim to increase this sample dKane et al.ll2009l) . 
the overall yield is expected to be small. 

Because RV surveys generically require spectroscopic ob- 
servations that are observationally expensive and must be 
obtained in series, it is more efficient to discover transiting 
planets around the much more abundant fainter stars by first 
searching for the photometric transit signal, and then follow- 
ing these up with targeted RV observations to eliminate false 
positives and measure the planet mass. However, in order to 
compensate for the rarity and low duty cycle, many stars must 
be monitored over a long time baseline. Photometric transit 
surveys that target brighter stars therefore require larger fields 
of view. Most of the original transit surveys had fields of view 
and exposure times that were optimized to detect planets or- 
biting stars with V > 10. Indeed, only ~ 20 transiting planets 
orbiting stars with V < 10 are currently known (^40 with 

V < 11). Of those with V < 10, ~ 40% were originally de- 
tected by RV surveys. 

The Kilodegree E xtremely Little Teles cope-North (KELT- 
North) transit survey (Pep per et al. 120071) was designed to de- 
tect giant, short-period transiting planets orbiting the bright- 
est stars tha t are not readily accessible to RV surveys. 
Peppe ret al. I d2003) determined the optimal hardware setup 
specifically to detect transiting planets orbiting stars with 

V ~ 8 - 10, and based on the specified design requirements in 
that paper, the KELT-North survey telescope system was con- 
structed using off-the-shelf, high-end consumer equipment. 
In fact, as the current detection demonstrates, KELT has ex- 
ceeded its design goals, and is sensitive to transiting systems 
in some favorable cases down to V ~ 12. 

In addition to the goal of filling in the magnitude gap be- 
tween radial velocity and other transit surveys, the KELT- 
North survey also has the potential to detect fainter systems 
with V > 10 that are in the magnitude range of previous sur- 
veys, but were missed or overlooked for various reasons. The 
detection discussed in this paper is an example of this oppor- 
tunity. Here the fact that the KELT-North survey is only now 
starting to vet candidates, more than eight years after the first 
candidates were announced by other transit surveys, can be 
seen an advantage. In particular, previous surveys have es- 
tablished the existence of massive brown dwarf companions 
dDeleuil et alj 120081: lirwin et al.1 1201 Ot iBouchv et all 1201 lah 



Johns on et al. 



2011b IBouchv et alJl201 lbl) . and have demon- 
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strated the feasibility of detecting low-mass compani ons to 
hot, rapidly rotating stars (ICollier Cameron et al.ll2010b . Par- 
tially in response to these results, the KELT-North survey de- 
liberately broadened our search targets to include hot and/or 
rapidly -rotating stars, which were previously neglected by 
many transit surveys. The evolving perception of what kinds 
of stars constitute viable transit search targets played an inter- 
esting role in the discovery of KELT- lb, as discussed in 33.21 

The KELT-North survey has been collecting data since 
September 2006, and has acquired a sufficient number of 
high-quality images to detect transit candidates. We have 
been systematically identifying and vetting transit candidates 
since February 201 1, and in this paper we report our first con- 
firmed low-mass transiting companion, which we designate 
KELT-lb. KELT-lb has a mass of ~ 27 Mj up , and we will 
therefore follow convention and refer to it as a 'brown dwarf 
throughout the majority of this paper. However, as we dis- 
cuss in 36.11 we are, in fact, agnostic about its true nature and 
therefore how it should be categorized. 

The outline of this paper is as follows. In order to intro- 
duce the survey and provide the appropriate context for our 
discovery, in $2] we summarize the properties of the KELT- 
North survey and our procedure for candidate selection. In 
$3] we review the observations of KELT-1, starting with the 
properties of the candidate in the KELT-North data, and then 
summarize the follow-up photometry, spectroscopy, and high- 
contrast imaging, ^describes our analysis and characteriza- 
tion of the host star and its substellar companion. In §|6] we 
provide a speculative discussion of the possible implications 
of this unique system for theories of the emplacement and 
tidal evolution of short-period substellar companions, mod- 
els of the structure and atmosphere of brown dwarfs, and the 
demographics of substellar companions to stars. We briefly 
summarize in Sj7] 

2. THE KELT-NORTH SURVEY 

Because this is the first paper from the KELT-North sur- 
vey, we describe the survey, selection criteria, and follow- 
up observations and reduction methodology in some detail. 
Readers who are not interested in these details, but are rather 
primarily interested in the properties and implications of the 
KELT-lb system, can skip to ^5] 

2.1. KELT-North Instrumentation and Survey Strategy 

The KELT-North survey instrument consists of a collec- 
tion of commerciall y-available equipme nt, chosen to meet 
the requirements of Pep per et al. I (120031) and tuned to find 
the few brightest stars with transiting planets in the Northern 
sky. The optical system consists of an Apogee AP16E (4K x 
4K 9/im pixels) thermo-electrically cooled CCD camera at- 
tached using a custom mounting plate to a Mamiya camera 
lens with a 80mm focal length and 42mm aperture (f/1.9). 
The resultant field of view of the detector is 26° x 26° at 
roughly 23"per pixel, allowing simultaneous observation of 
nearly 40,000 stars in typical high Galactic latitude fields. 
The medium-format image size is markedly larger than the 
CCD detector (which measures 37 x 37mm) which greatly re- 
duces the severity of vignetting across the large field of view. 
At the same time, the small aperture permits longer expo- 
sures, which improve observing efficiency (assuming fixed 
camera read-out time). A Kodak Wratten #8 red-pass filter 
is mounted in front of the lens to further reduce the impact 
of atmospheric reddening (which primarily affects blue wave- 
lengths) on our photometry. The resultant bandpass resembles 



a widened Johnson-Cousins R-band. This optical system is 
mounted atop a Paramount ME robotic mount from Software 
Bisque on a fixed pier at Winer Observatory in Sonoita, AZ 
(Latitude 31° 39'56.08"N Longitude 110° 36' 06.42"W, ele- 
vation 1515.7 meters). See Pepper et al. (2007) for additional 
details about the system hardware. 

The primary KELT-North transit survey consists of 13 fields 
centered at 3 1 .7 declination, spanning all 24 hours of right as- 
cension. Including the slight overlap between fields, the total 
survey area is « 40% of the Northern sky. Survey observa- 
tions consist of 150-second exposures with a typical per- field 
cadence of 15-30 minutes. The KELT-North telescope has 
been collecting survey data in this manner since September 
2006, and to date has acquired between 5000 and 9300 images 
per field. Given this quantity of data and the typical achieved 
photometric precision of ~ 1% for V < 11, the KELT-North 
survey is able to detect short-period giant transiting planets or- 
biting most FGK stars with magnitudes from saturation near 
V ~ 8 down to V ~ 12. 

2.2. KELT-North Pipeline 

Relative photometry is generated from fl at-fielded images 
using the ISIS im age subt raction package dAlard & L upton 
ll998tlAlardll2000l see also lHartman et al.ll2004l) . in combina- 
tion with point -spread function fitti ng using the stand-alone 
DAOPHOT II ( IStetsod[l98ll[T990h . Although highly effec- 
tive, the image subtraction procedures are highly computer- 
intensive. To improve reduction performance, the default ISIS 
scripts were modified to facilitate distributed image reduction 
across many computers in parallel. ISIS operation in this fash- 
ion permits thorough exploration of various reduction param- 
eters, which would be intractable if executed serially. Other 
elements of the ISIS reduction package have also been mod- 
ified or replaced with more robust alternatives. For example, 
the standard ISIS source-identification routines and utilities 
are ill-equipped to deal with the nature and ubiquity of aber- 
rations in KELT-North images. In response, we have replaced 
the ISIS 'extract' utility w ith the popular SExtractor program 
dBertin & Arnoutsl[T9 96). A more complete explanation of 
these modifications and driver scripts that implement them are 
available online 17 . 

2.3. KELT-North Candidate Selection 

Once we have the light curves created by ISIS for all of the 
DAOPHOT-identified point sources in the reference image, 
we begin a series of post-processing steps before doing the 
initial candidate selection. To begin, we convert the ISIS light 
curves from differential flux to instrumental magnitude using 
the results of the DAOPHOT photometry on the reference 
image. We also apply 5a iterative clipping to all of the light 
curves at this stage; this typically removes ~ 0.6% of the data 
points. All of the uncertainties for the converted and clipped 
light curves in a given field are then scaled as an ensemble. 
The scaling is chosen such that the x 2 /dof = 1 for the main 
locus of the light curves on a magnitude versus x 2 /dof plot. 
Typically this scaling is around a factor of 1 .2, implying that 
the uncertainties are somewhat underestimated. 

We next attempt to match all of the DAOPHOT-identified 
point sources in the reference image to stars in the Tycho-2 
catalog. We obtain a full-frame WCS with sub- pixel accu- 
racy on our reference frame using Astrometry.net (lLang et af] 

17 http://astro.phy.vanderbilt.edu/~siverdrj/soft/is3/index.html 
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2010). Using this solution, we match stars by taking the 
closest Tycho-2 entry within 45". This typically generates 
matches for 98% of the Tycho-2 stars within each field. A 
successful Tycho-2 match also will provide a 2MASS ID. We 
use the proper motions and JHK apparent magnitudes from 
these two catalogs. 

With this catalog information, we next identify and exclude 
giant stars by means of a red uced proper motion (Hj) dia- 
gram (iGould & Morgan 2003); F ollowing the specific pre- 
scription of lColl ier Came ron et al.l (l2007bl) . we place each of 
our matched stars on a J vs. Hj plot. We compute the reduced 
proper motion of a star as 



Hj = 7+51og(/i/mas yr ) 



(1) 



and determine the star to be a giant if it fails to satisfy 

Z7, > -141 .25(7- //) 4 +473 . 1 8(7- i/) 3 

-583 .6(7- Hf + 313 .42(7- H) - 43 .0 (2) 



This process leaves us with anywhere from 10,000 to 30,000 
catalog-matched putative dwarf stars and subgiants (hereafter 
dwarfs) per field, depending primarily on the location of the 
field relative to the Galactic plane. 

The dwar fs are then run throu gh the Trend Filtering Algo- 
rithm (TFA. lKovacs et ai1l2005l) 18 to reduce systematic noise. 
We select a new set of detrending stars for each light curve by 
taking the 150 closest stars - outside of a 20 pixel exclusion 
zone centered on the star being detrended - that are within 
two instrumental magnitudes of the star being detrended. 

KELT'S Paramount ME is a German Equatorial mount, 
which requires a "flip" as it tracks stars past the meridian. 
Therefore, the optics and detector are rotated 180 degrees 
with respect to the stars between observations in the Eastern 
and Western hemispheres, and detector defects, optical dis- 
tortions, PSF shape, flat fielding errors, etc., for a given star 
can be completely different. This requires us to treat obser- 
vations in the East and West essentially as two separate in- 
struments. Thus the preceding steps (magnitude conversion, 
error scaling, dwarf identification, TFA) are each performed 
separately on the East and West images of each field. After 
the dwarf stars in the East and West have been run through 
TFA, we then combine the two light curves of each target into 
one East+West light curve. We first match stars from the East 
and the West pointings by their Tycho IDs, and then deter- 
mine the error- weighted scaling factor of the Western light 
curve needed to match the error-weighted mean instrumental 
magnitude of the East light curve. 

All of the light curves from the matched Tycho dwarf stars 
in a field are given an internal ID. We next search the com- 
bined East+West light curves of the dwarfs for transit-like 
signals using the b ox-fitting least squares algorithm (BLS; 
iKovacs et alJl2002T) . We use a custom version of BLS modi- 
fied to skip over integer and half integer multiples of the side- 
real day to reduce the effect of spurious signals due to diurnal 
systematics and their aliases on the BLS statistics. We per- 
form selection cuts along six of the statistics that are output 
by the VARTOOLS implementation of the BLS algorithm: sig- 
nal detection efficiency SDE, signal to pink noise SPN 19 , the 
fraction of transit points from one night f\ n , depth 5, the ra- 
tio of A\ 2 for the best transit model to best inverse transit 

18 We used the v ersions of TFA and BLS (described later) found in the 
VARTOOLS package <Hart man et al. 2 003) . 

" See Kovacs et al. 12002) and Hartman et al. 1 2009), respectively, for the 
definitions of SDE and SPN 



TABLE 1 

KELT-NORTH BLS CANDIDATE SELECTION CRITERIA 



Signal Detection Efficiency SDE>7.0 
Signal to Pink-noise SPN>7.0 
Fraction from one night f\ n < 0.8 



Depth 
X 2 ratio 



6 < 0.05 



Duty cycle q < 0. 1 



> 1.5 



model Ax 2 /Ax- dBurke et al.ll2006l) . and the fraction of the 
orbit spent in transit or duty cycle q. In order to determine the 
appropriate threshold values for these statistics, we injected 
realistic transit signals with a range of properties into a large 
sample of light curves, and then attempted to recover these 
using the BLS algorithm. We then determined the values of 
these statistics that roughly maximize the overall detection ef- 
ficiency while minimizing the number of spurious detections. 
The final adopted values are given in Table [TJ 

In addition to the cuts we make on the BLS statistics, we 
also impose restrictions on the effective temperature and in- 
ferred density of the candidate host stars. For the temperature, 
we require that T e g < 7500K. We calculate the stellar effective 
temperature of each candidate from it s 2MASS J-K colors . 
We used the Yonsei-Yale isochrones (Dem arque et al.ll2004l) 
at 5 Gyr with solar metallicity and no alpha enhancement to 
create a simple polynomial fit for T s s as a function of J-K: 

log T eff = 3. 94808 - 0.7353(7- tf) 

+1.0116(7-/C) 2 -0.8334(7-/i:) 3 . (3) 

As we have conducted our follow-up spectroscopy, we have 
found that this relation generally predicts T e g to within ~ 
100K for T eB < 7000K and to within - 300K for stars with 
r eff = 7000-7500K. 

We also require that the stellar density, p„, as inferred from 
the BLS transit fit to the KELT-North light curve, to be within 
1.0 dex of the stellar density calculated for each star using its 
J-K colors, assuming the star is on the main sequence. A 
large disparity in the observed versus the calculated density 
is indicative of a blend or o f a giant that made it through th e 
reduced proper motion cuts (ISeager & Mallen-Orn elas 2003 ). 
Again using the Yonsei-Yale isochrones at 5 Gyr with solar 
metallicity and no alpha enhancement, we made a fit for den- 
sity as a function of J-K: 

log(p*,c a ic/po) = -l-00972 + 2.82824(/- J fT) 

-1.19772(7-/C) 2 . (4) 

We require that this value be within 1.0 dex of the stellar den- 
sity we calculate from the KELT-North lightcurve 



logp >t , obs = log 



Gir 2 q 3 P 2 



(5) 



where P and q are the orbital period and duty cycle (tran- 
sit duration relative to the period) as returned by BLS. This 
equation assumes circular orbits and that the companion mass 
is much smaller than the host star mass, Mp <C M* . Also, 
because BLS does not attempt to fit for the ingress/egress 
duration, and furthermore KELT-North data typically do not 
resolve the ingress or egress well, we are not able to deter- 
mine the transit impact parameter and thus the true stellar ra- 
dius crossing time. Equation [5] therefore implicitly assumes 
an equatorial transit, and so formally provides only an upper 
limit to the true stellar density. For a transit with an impact pa- 
rameter of b = 0.7, the true density is ~ 0.5 dex smaller than 
that inferred from Equation 
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All of the light curves that pass these selection crite- 
ria are designated as candidates, and a number of addi- 
tional diagnostic tests are then performed on them , includ- 
ing Lomb-Scarele (LS.lLombl|1976l:IScarglelll982h and AoV 
(ISchwarzenberg-Czernyl 119891 iDevorl 120051) periodograms. 
The results of these tests, the values of the BLS statistics, the 
light curves themselves, as well as a host of additional infor- 
mation, are all collected into a webpage for each candidate. 
Members of the team can then use this information to vote 
on the true nature of the candidate (probable planet, eclips- 
ing binary, sinusoidal variable, spurious detection, blend or 
other). All candidates with at least one vote for being a prob- 
able planet are then discussed, and the most promising are 
then passed along for follow-up photometry, reconnaissance 
spectroscopy, or both. 

3. OBSERVATIONS 

3.1. KELT-North Photometry, Candidate Identification, and 
Vetting Overview 

KC20C05168 emerged as a strong candidate from the anal- 
ysis of the combined light curves from stars in the overlap 
area between fields 1 and 13. The KC20C05168 light curve 
contains 8185 epochs distributed over ~ 4.2 years, between 
UT October 25, 2006 and UT December 28, 2010, with a 
weighted RMS of 9.8 millimagnitudes (mmag). This RMS is 
typical for KELT-North light curves of stars with this magni- 
tude (V ~ 10.7). A strong BLS signal was found at a period of 
P ~ 1 .2175 days, with a depth of d ~ 3.8 mmag, and detection 
statistics SPN=8.53, SDE= 12.41, q = 0.09, A X 2 /A X 2 - = 2.06, 
and log(p* ,obs/p* ,cai) = -0.06. The phased KELT-North light 
curve is shown in Figure [T] A significant signal also appeared 
in SuperWASP data (B utters et al.l 1201 Oft of this star at the 
same period. The KELT-North data exhibit some evidence 
for out-of-transit variability at the mmag level and exhibit 
some relatively weak peaks in the LS and AoV periodograms, 
but we did not consider these signals to be strong enough 
to warrant rejection of the candidate. In addition, the depth 
of the photometric transit signal in the original KELT-North 
light curve is substantially smaller than we find in the high- 
precision follow-up data (see 93.4b . Further analysis indicates 
that the out-of-transit variability and smaller depth were likely 
due to a minor problem with the original data 

Based on the strength of the K20C05168 signal, the esti- 
mated effective temperature of the host star of T e e ~ 6500K, 
and the fact that the star was sufficiently isolated in a DSS 
image, we submitted the candidate for reconnaissance spec- 
troscop y with the Till inghast Reflector Echelle Spectrograph 
(TRES: iFCriszl 12008ft on the 1.5m Tillinghast Reflector at 
the Fred Lawrence Whipple Observatory (FLWO) on Mount 
Hopkins in Arizona. The first observation on UT November 
9, 2011 at the predicted quadrature confirmed the T e ff esti- 
mate of the star, and also demonstrated that it was a slightly 
evolved dwarf with logg ~ 4, and that it was rapidly rotating 
with vsin/* ~ 55 km s" 1 . A second observation was obtained 
on UT November 1 1, 201 1 separated by ~ 1 .9 days, or ~ 1 .54 
in phase, from the first observation and thus sampled the op- 
posite quadrature. The two observations exhibited a large and 
significant radial velocity shift of ~ 8 km s" 1 , consistent with 
a brown dwarf companion. 

Efforts to obtain photometric follow-up during the primary 
transit and secondary eclipse were then initiated. Concur- 
rently, additional spectra with TRES were taken to character- 
ize the spectroscopic orbit. In addition, we obtained adaptive 
optics imaging of the target to search for close companions. 
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FIG. 1.— The KELT-North light curve of KELT-1 phased to the BLS 
determined period of P = 1.2175 days is shown in the grey points. The black 
points show the data binned 0.02 in phase. 

Finally, once we were fairly confident that the signals were 
due to a low-mass transiting companion, we obtained contin- 
uous spectroscopic time series with TRES during the primary 
transits on UT December 21, 2011 and UT January 7, 2012 
for the purposes of measuring the Rossiter-McLaughlin (RM) 
effect. All of these observations are described in greater detail 
in the subsequent sections and summarized in Table [2] 

3.2. Previous Identification of the Photometric Candidate by 

HATNet 

KELT- lb was also recognized as a photometric transiting 
planet candidate by the HATNet project, based on observa- 
tions obtained in 2006. In September 2009 the candidate 
HTR162-002 was forwarded to D. Latham's team for spectro- 
scopic follow up. An initial observation with TRES confirmed 
that the target was a late F main-sequence star, as expected 
from 2MASS color 7-^=0.245. The synthetic template with 
r e ff = 6250K and logg = 4.0 and assumed solar metallicity, 
gave the best match to the observed spectrum. However, that 
first TRES spectrum also revealed that the star was rotating 
rapidly, with vsin/* = 55 km s" 1 . At that time, D. Latham's 
team routinely put aside candidates rotating more rapidly than 
about vsin/* = 30 km s , arguing that it would not be pos- 
sible to determine velocities with a precision sufficient for an 
orbital solution for a planetary companion. 

HTR162-002 remained on the HATNet "don't observe 
with TRES" list until it was independently rediscovered 
by the KELT-North team and was forwarded as candidate 
KC20C05168 to D. Latham's team in November 2011 for 
spectroscopic follow up with TRES. During the intervening 
26 months, there were two relevant developments in the pro- 
cedures and tools used by Latham's team, both resulting from 
contributions by L. Buchhave. The first development, enabled 
by convenient tools in the observing website, was the prac- 
tice of observing new candidates only near opposite quadra- 
tures, according to the discovery ephemeris and assuming cir- 
cular orbits. The second development was a much improved 
procedure for deriving radial velocities for rapidly-rotating 
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stars, initially motivated by the Kepler dis covery of hot whit e 
dwarfs transiting rapidly-rotating A stars dRowe et al.ll2010h . 
As it turned out, the second observation of KC20C05 168 with 
TRES described above was taken before the first observation 
was reviewed, so the candidate was not relegated to the re- 
jected list due to its rapid rotation before the opposite quadra- 
ture was observed. When the results were reviewed after the 
second observation, the evidence for a significant radial ve- 
locity shift between the two observations was obvious, despite 
the rapid rotation, therefore suggesting that the unseen com- 
panion was probably a brown dwarf, if not a giant planet. 

It should also be recognized that over the 26 months since 
the first observation of HTR 162-002, the attitude against pur- 
suing rapidly rotating stars as possible hosts for transiting 
planets had gradually softened among the exoplanet com- 
munity. An important development was the demonstration 
that slowly-rotating subgiants that have evolved from rapidly- 
rotating main-sequence A stars do occasionally show the 
radial -velocity signatures of orbiting planetary companion 
(e.g.. iJohnson et aii r2007). A second insight came from the 
demonstration that the companion that transits the rapidly - 
rotating A star WASP-33 must be a planet, using Doppler 
imaging ( Collier Ca meron et al.ll2010h . Finally, the discovery 
of transiting brown dwarf companions suggested the possibil- 
ity of detecting their large amplitude RV signals even when 
they orbit stars with large vsin/* and thus poor RV precision. 

In the early days of wide-angle photometric surveys for 
transiting planets, Latham's team had established procedures 
for handling candidates forwarded for spectroscopic follow up 
by more than one team. Such duplications were fairly com- 
mon, and the goal was to assign credit to the initial discov- 
ery team, which was especially important in an era when few 
transiting planets had been confirmed. By the time it was no- 
ticed in mid December 201 1 that KC20C05168 was the same 
as HTR 162-002, the KELT-North team already had in hand a 
convincing orbital solution from TRES and high-quality light 
curves from several sources, confirming that KELT- lb was in- 
deed a substellar companion. 

3.3. Spectroscopy from FLWO/TRES 

A total of 8 1 spectra of KELT- 1 were taken using the TRES 
spectrograph on the 1.5m Tillinghast Reflector at FLWO. 
These were used to determine the Keplerian parameters of 
the spectroscopic orbit, measure bisector variations in order 
to exclude false positive scenarios, measure the spectroscopic 
parameters of the primary, and measure anomalous RV shift 
of the stellar spectral lines as the companion transits i n front 
of the rapidly-rotating ho st star, i.e., the RM effect dRossiterl 
1 1924 lMcLaughlinlll924l) . The TRES spectrograph provides 
high resolution, fiber- fed echelle s pectroscopy over a band- 
pass of 3900-8900 dFureszl 120081) . The observations ob- 
tained here employed the medium fiber for a resolution of 
R ~ 44,000. The dat a were reduced and an alyzed using the 
metho ds described in lOuinn et al.l d2012l) and lBuchhave et all 

d2oToh . 

A subset of six spectra were combined in order to determine 
the spectroscopic parameters of the host star using the Spec- 
tral Parameter Classification (SPC) fitting program (Buchhave 
etal., in preparation). SPC cross-correlat es the observed spec- 
trum against a grid of synthetic Kurucz (iKuruczll 1979b spec- 
tra. This analysis yielded T eS = 6512 ± 50K, logg = 4.20 ± 
0.10, [Fe/H]= 0.06 ± 0.08, and vsin/* = 55.2 ± 2km s" 1 . 
These parameters were used as priors for the joint global fit 
to the RV, RM, and photometric data as described in §15.21 




860 880 900 920 940 960 980 
BJD TnR - 2455000 



FIG. 2. — (Top panel) The black points with uncertainties show the mea- 
sured RVs for KELT-1 as a function of time in BJDtdb- The barycentric 
veloc ity of the system, as determined from the model fit shown in red (see 
j|5. 2t . has been been subtracted from the data. (Bottom panel) The residuals 
from the model fit. 

Spectra were taken at a range of phases in order to char- 
acterize the spectroscopic orbit and search for bisector span 
variations indicative of a blend. One of these spectra hap- 
pened to be taken during a primary transit on UT 201 1-1 1-18, 
and so was not used in the analysis because it is likely to be 
affected by the RM effect. The RV and bisector data for the 
remaining 23 spectra are listed in Table [7] These observations 
span ~ 88 days from UT 2011-11-09 through UT 2012-02- 
05. The uncertainties on the listed radial velocities have been 
scaled by a factor of 1 .214 based on an independent fit to these 
data, as described in §15.21 The scaled median RV uncertainty 
is ~ 230 m s" 1 . The uncertainties in the bisector measure- 
ments have not been scaled. The median bisector uncertainty 
is ~ 110m s" 1 . 

Time series spectroscopy was obtained with TRES on two 
different nights of primary transits in order to measure the 
spin-orbit alignment of the companion via the RM effect. Fif- 
teen observations were obtained on UT 201 1-12-21 and forty- 
two observations on UT 2012-01-07. Conditions were rela- 
tively poor for the first run, resulting in a factor ~ 2 larger un- 
certainties and incomplete coverage of the transit. We there- 
fore decided not to include these data in our final analysis, 
although we confirmed that this has no effect on our final in- 
ferred parameters. The RV and bisector data for the RM run 
on UT 2012-01-07 are listed in Table The RV uncertainties 
have been scaled by a factor of 0.731, also based on the global 
model fit described in §15.21 We note that the majority of the 
X 2 for these data are due to a few outliers. The median scaled 
RV uncertainty is ~ 160 ms" 1 . The bisector uncertainties 
were not scaled. 

All of the RV and bisector measurements used in the sub- 
sequent analysis are shown as a function of epoch of observa- 
tion in BJDtdb in Figure [3] The measurements phased to the 
best-fit companion period from the joint fit to photometric and 
RV data are shown in Figure [3] demonstrating the very high 
signal-to-noise ratio with which the RV signal of the compan- 
ion was detected, and the good phase coverage of the orbit. 
A detail of the RV data near the time primary transit with the 
orbital (Doppler) RV signal removed is shown in Figure |4] 
showing the clear detection of the RM effect and a suggestion 
that the orbit normal is well aligned with the projected stellar 
spin axis. 

Finally, we determined the absolute radial velocity of the 
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FIG. 3. — The black points with uncertainties show the measured RVs 
for KELT- 1 relative to the barycentric velocity of the system, pha sed t o the 
best-fit period as determined from the model fit shown in red (see j|5.2t . The 
phases are normally referenced to the time of periastron (7>), but have been 
shifted such that a phase of 0.25 corresponds to the time of inferior conjunc- 
tion Tc or primary transit. RV data near this phase show deviations from 
the Keplerian expectation due to the RM effect, which was included in the 
model. (Middle panel) The residuals of the RV data from the model fit. (Bot- 
tom panel) Bisector spans as a function of phase. 
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FIG. 4. — (Top panel) The black points with uncertainties show the mea- 
sured RVs relative to the barycentric velocity of the system for KELT- 1 , as a 
function of the time since primary transit Tc, for data taken near Tc- The Ke- 
plerian RV variation as determined from the best-fit model has been removed 
from both the data and model. Data taken within ~ 1.4 hours of Tc occur 
during primary transit, and are thus strongly affected by the RM effect. The 
shape of the RM signal indicates that the projected obliquity of the host star 
with respect to the orbit is small. (Bottom panel) The residuals of the data to 
the RM fit. 

system barycenter using a simple circular orbit fit to radial 
velocities determined from the full set of spectra, which were 
determined using a single order near the Mg b line. (Note 
that the relative RVs used for determining the orbit were de- 
termined using the full, multi-order analysis of the spectra.) 
The zero point correction to these velocities were determined 
using contemporaneous monitoring of five RV standard stars. 
The final value we obtain is 7„b s = -14.2 ± 0.2 km s" 1 , where 
the uncertainty is dominated by the systematic uncertainties in 
the absolute velocities of the standard stars. This zero point, 
along with the global fit to the data in 35.21 were used to place 
the instrumental relative radial velocities on an absolute scale. 
Therefore, the RVs listed in Tables [7] and [8] are on an absolute 
scale. 

3.4. Follow-Up Photometry 



1.20 



c 
to 

CO 

c 
o 

o 



X 

_2 

X> 
CD 
N 

E 



1.15 



1.10 



1.05 



1.00 



mi in in iiii|iii inii|ii in i|i iiniii| in mi in ii 

PvdKO UT 2011-12-03 



ULMOUT 201 1-12-03 



HAO UT 2011-1 2-10 



FLWOUT 201 1-12-1 6 t 



ULMOUT 2011-12-31 



FLWOUT 2012-01-07 




I 



-3 



-10 12 
Time - T c (hrs) 



FIG. 5. — The black points show the relative flux as a function of time 
from primary transit ( Tc) for the six sets of follow-up observations of primary 
transits we analyze here. The data sets are labeled and summarized in Table 
f2] The data are normalized by the fitt ed o ut-of-transit flux, and a linear trend 
with airmass has been removed (see jj5.2t . In addition, an arbitrary offset has 
been applied to each light curve for clarity. For each observation, we plot the 
data above and the resid uals below. In all cases, the red lines show the model 
fit from the analysis in jj5.2| 

We obtained high-precision follow-up photometry of 
KELT-1 in order to confirm the K20C05168 transit signal, 
search for evidence of a strongly wavelength-dependent tran- 
sit depth indicative of a stellar blend, and search for evidence 
of a secondary eclipse. Also, these data enable precision mea- 
surements of the transit depth, ingress/egress duration, and 
total duration, in order to determine the detailed parameters 
of the KELT-1 system. In all, we obtained coverage of 9 
complete and 4 partial primary transits, and two complete and 
one partial secondary eclipse, using six different telescopes in 
all. Many of these data were taken under relatively poor con- 
ditions and/or suffer from strong systematics. We therefore 
chose to include only a subset for the final analysis, including 
six of the primary transits and the three secondary eclipses. 
In the following subsections, we detail the observatories and 
data reduction methods used to obtain these data. The dates, 
observatories, and filters for these data sets are summarized 
in Table [2] The light curves for the primary transits are dis- 
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FIG. 6. — The grey points show the combined i and Pan-STARRS-Z rela- 
tive photometry of KELT-1 as a function of the predicted time of secondary 
eclipse of KELT-lb (7^), obtained from the observatories listed in Table|2] 
The data have been corrected for a lin ear t rend with airmass and normalized 
by the zero point of the linear fit (see jj5.4l . The larger circles with error bars 
show the binned observations. Note we do not fit to the binned data; these are 
shown for the purposes of illustration only. The over plotted example light 
curve is the secondary eclipse depth we would expect if KELT-lb had a ge- 
ometric albedo of A g = 0. 1 and instantaneously reradiated its incident stellar 
flux (/' = 2/3). We would have detected this event with a confidence level of 
> 95%. 
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FIG. 7. — (Top panel) The black points show the six data sets displayed 
in Figure[5] combined and binned in 5 minute intervals. Since these data sets 
were taken with different filters and have different systematics, we do not use 
this combined light curve for analysis, but rather show it for the purposes of 
illustrating the overall quality and statistical constraining power of the en- 
semble of follow-up light curves. The red curve shows the six transit models 
for each of the individual fits combined and binned in 5 minute intervals the 
same way as the data. (Bottom panel) The residuals of the binned light curve 
from the binned model in the top panel. 

played in Figure [5]and the data are listed in Tables [9]through 
n~3l whereas the light curves for the secondary eclipse are dis- 
played in Figure[6]and the data are listed in Tables [T5lthrough 

El 

3.4.1. Peter van de Kamp Observatory (PvdKO) 

Data on the primary transit starting on UT 201 1-12-03 were 
acquired with the 0.6-meter, f/7.8 Ritchey-Chretien telescope 
at the Peter van de Kamp Observatory at Swarthmore College 



(Swarthmore, PA, USA). The telescope is equipped with an 
Apogee U16M CCD with 4096x4096 9-micron pixels, giv- 
ing a field of view 26 arcminutes on a side. Available fil- 
ters are 50mm-square Johnson-Cousins UBVR C I C and SDSS 
ugriz, both from Omega Optical. The telescope is autogu- 
ided to minimize photometric errors from imperfect flatfield- 
ing, keeping the centroid shift of each star to within typically 
3-4 pixels over the course of a night. The observations used 
here were obtained with the i filter and used 2x2 binning, giv- 
ing a binned pixel scale of 0.76"/ pixel. 

The data were reduced in IRAF using standard procedures 
for flat-fielding (with twilight sky flats) and dark and bias sub- 
traction. Aperture photometry was performed, and then dif- 
ferential magnitudes for the target star were calculated using 
an ensemble of comparison stars in the same field, chosen to 
minimize the scatter in the final light curve. 

3.4.2. University of Louisville Moore Observatory (ULMO) 

Data on the primary transits starting UT 2011-12-03 and 
2011-12-31, and on the secondary eclipses starting 2011-12- 
02 and 2012-01-04, were obtained with the University of 
Louisville Moore Observatory RC24 telescope (MORC24) 
located near Brownsboro, Kentucky. MORC24 is a RC Op- 
tical Systems Ritchey-Chretien 0.6 m telescope on an equa- 
torial fork mount. The telescope is equipped with an Apogee 
U16M 4096x4096 pixel CCD camera which has a focal plane 
scale of 0."39 pixel" 1 and a field of view (FOV) of 26 .'3 x26'.3. 
The UT 2011-12-03 and 2011-12-31 data were obtained us- 
ing an Astrodon Photometries Sloan r filter, while the other 
two sets of data were obtained using an Astrodon Photomet- 
ries Sloan i filter. The MORC24 mount has excellent free- 
running tracking, so we did not use a separate guide camera. 
Instead, minor telescope pointing corrections are made after 
each exposure by comparing the CCD pixel coordinates of the 
centroid of the target star to its initial position on the CCD. 
KELT-lb was held to within 3-4 pixels of the starting position 
on the CCD throughout each observing session. Since KELT- 
lb is separated from its nearest detectable neighbor in DSS2 
imagery by ~ 18", we were able to defocus the telescope to 
allow for longer exposures without the risk of blending from 
the neighbor star. An exposure time of 100 s was used for 
all observations, resulting in a 120 s cadence when combined 
with the 20 s CCD readout time. 

We used AstroImageJ (Collins & Kielkopf 2012, in prepa- 
ration) to calibrate the image data. The algorithm includes 
bias subtraction, CCD non-linearity correction, dark subtrac- 
tion, and flat-field division. AstroImageJ was also used to per- 
form aperture photometry using a circular aperture. An aper- 
ture size and an ensemble of comparison stars in the same 
field were chosen to minimize the scatter in the final light 
curves. AstroImageJ provides the option to use a standard 
fixed radius aperture or a variable radius aperture based on 
the measured FWHM of the target star in each image of the 
series. When a star is well separated from other stars, the 
variable aperture option tends to reduce photometric scatter 
under observing conditions that result in significant changes 
to the PSF during the observing session. The variable aper- 
ture produced optimal results for all four MORC24 KELT-lb 
light curves. 

For the observations starting on UT 2011-12-02, cirrus 
clouds were present during the first half of the observations, 
and airmass ranged from 1.16 at the start of observations to 
3. 19 at the end. For the observations starting on UT 201 1-12- 
04, skies were clear until clouds moved in about 30 minutes 
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after ingress. The clouds cleared just prior to egress, how- 
ever, sky transparency remained highly variable until about 
an hour after egress. Airmass ranged from 1.05 at the be- 
ginning of observations to 1.40 at the end. Although guid- 
ing was maintained through the cloud cover, data during that 
time have been removed. For the observations starting on UT 
20 11-12-31, skies were clear with minimal variations in trans- 
parency. Airmass ranged from 1 .00 at the beginning of obser- 
vations to 2.17 at the end. For the observations on UT 2012- 
01-04, cirrus clouds were present during the second half of 
the observations, and airmass ranged from 1.03 at the start of 
observations to 1 .96 at the end. 

3.4.3. Hereford Arizona Observatory (HAO) 

Data on the primary transit starting UT 2011-12-10 were 
obtained at the Hereford Arizona Observatory, HAO (obser- 
vatory code G95 in the IAU Minor Planet Center). This is 
a private observatory in Southern Arizona consisting of a 14- 
inch Meade LX-200 GPS telescope equipped with a SBIG ST- 
10XME CCD, a focal reducer and a 10-position filter wheel 
with SDSS filters ugriz. The telescope and dome are oper- 
ated via buried cables, permitting automation of observing 
sessions. Calibrations usually employ a master fiat frame ob- 
tained during dusk prior to the observing session. The field- 
of-view (27 x 18 arcminutes) is sufficient for the use of ap- 
proximately two dozen stars as candidates for reference in a 
typical field. The observations reported here were obtained 
with the i filter. 

The data were reduced and a light curve was generated as 
follows. An artificial star was inserted in each image before 
photometry readings for the purpose of monitoring smooth 
extinction as well as extra extinction events caused by thin 
clouds, dew formation, and atmospheric seeing degradations 
that could swell the PSF beyond the photometry aperture cir- 
cle. Photometry magnitude readings were made by the pro- 
gram Maxim DL and imported to a spreadsheet, where sev- 
eral steps of manual reduction were completed. The first was 
to solve for an extinction model (including a temporal extinc- 
tion trend) based on the sum of all candidate reference star 
fluxes versus air mass. Second, subsets of reference stars were 
evaluated for suitability, by toggling individual stars "on and 
off" in order to determine the subset that minimize the RMS 
scatter in the target star light curve. 

Finally, the light curve for the target was fitted using a 
model for systematic effects and a transit signature. System- 
atics were represented by a temporal trend and air mass cur- 
vature (AMC). The AMC is caused by the target star having 
a color that differs from the flux- weighted color of the refer- 
ence stars. The transit parameters were depth, total duration, 
ingress/egress duration, and a parameter related to the stellar 
limb darkening. The solution was obtained by minimizing the 
X 2 of the fit. Outliers were identified using an objective re- 
jection criterion based on deviations from the model solution. 
Finally, the light curve is corrected for extinction and system- 
atic effects and scaled to the out-of-transit model flux. 

3.4.4. FLWO/KeplerCam 

Data on the primary transits on UT 2011-12-16 and 2012- 
01-07 were obtained with KeplerCam on the 1.2m telescope 
at FLWO KeplerCam has a single 4K x 4K Fairchild CCD 
with a pixel scale of 0.366 arcseconds per pixel, for a total 
FOV of 23.1 x 23.1 arcminutes. A full transit was observed 
on UT 2011-12-16 with clear conditions. Observations were 



obtained in the SDSS z filter with 30-second exposures. We 
also obtained a full transit on UT 2012-01-07 and observa- 
tions were obtained with the SDSS i filter with 15-second ex- 
posures. Clouds came in at the end of the transit and as a result 
there is some increased scatter in the out-of-transit baseline. 
The data w ere reduced using a light curve reduction pipeline 
outlined in Car ter et ail (1201 ll) which uses standard IDL tech- 
niques. 

3.4.5. Las Cumbres Observatory Global Telescope Network 

(LCOGTj 

Data on the secondary eclipse on UT 201 1-12-30 were ob- 
tained with the 2.0m Faulkes Telescope North (FTN) tele- 
scope, which is located on Haleakala on the island of Maui 
in Hawaii. The FTN telescope is part of the Las Cumbres Ob- 
servatory Global Telescope Network 20 . These observations 
were made using the 4K x 4K Spectral camera (Fairchild 
Imaging CCD486 BI) in bin 2x2 mode for a faster readout 
together with the PanSTARRS-Z filter. As scintillation noise 
becomes significant (>1 millimag) in exposures shorter than 
^30 sec for telescopes of this aperture, the exposure time was 
kept to 60 sec and the telescope defocused to avoid satura- 
tion of the target while ensuring sufficient signal-to-noise ra- 
tio in the comparison stars. These data were debiased, dark- 
subtracted and flat fielded by the LCOGT offline pipeline (de- 
veloped by the Astrophysics Research Institute at Liverpool 
John Moores) and aperture ph otometry was carrie d out using 
the stand-alone DAOPHOT II JStetsorjl987l[l99Ch . Differen- 
tial photometry was then computed using an ensemble of 15 
comparison stars. 

3.5. Keck Adaptive Optics Imaging 

To further assess the multiplicity of KELT-1, we acquired 
adaptive optics images using NIRC2 (PI: Keith Matthews) at 
Keck on UT 2012-01-07. Our observations consist of dithered 
frames taken with the K' (A c = 2.12^m) and H (A c = 1.65/im) 
filters. We used the narrow camera setting to provide fine 
spatial sampling of the stellar point-spread function. The total 
on-source integration time was 81 seconds in each bandpass. 

Images were processed by replacing hot pixel values, flat- 
fielding, and subtracting thermal background noise. No com- 
panions were identified in individual raw frames during the 
observations; however, upon stacking the images we noticed a 
point source (8<r) to the south-east of KELT- 1 . Figure[8]shows 
the final processed K' image. Inspection of the companion lo- 
cation showed that its separation from the star does not change 
with wavelength, demonstrating that it is not a speckle. This 
object is too faint and close to the primary to be detected with 
seeing-limited images. 

We performed aperture photometry to estimate the relative 
brightness of the candidate tertiary, finding AH = 5.90±0.10 
and AK' = 5.59 ± 0.12. An H - K' = 0.4 ± 0.2 color is 
consistent with spectral-t ypes M1-L0 (Le ggett et al.l 120021 
iKraus & "Hi llenbrand 2007). If the candidate is bound to 
KELT-1 and thus at the same distance of 262 ± 14pc and suf- 
fers the same extinction of Ay = 0.18 ±0.10 (see 35.11 ), then 
we estimate its absolute H magnitude to be Mh = 8.31 ±0.15, 
corresponding to a M4-5 spectral type, consistent with its 
color (see, e.g., the compilation of Kirkpatrick et al. 2012). 

We also measured an accurate position of the companion 
relative to the star by fitting a Gaussian model to each of the 

20 http://lcogt.net 
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FIG. 8.— A Keck AO image of KELT- 1 taken with NIRC2 on UT 2012- 
01-07 in the K' filter. North is up and East is to the left. A 0.5"bar is shown 
for scale. A faint companion with AK' = 5.59 ± 0. 12 located ~ 558 zt 1 mas 
to the southeast is clearly visible. 

point-spread function cores. After correcting for distortion in 
the NIRC2 focal plane 21 , and adopting a plate scale value of 
9.963 ± 0.006 mas pix -1 a nd instrume nt orientation relative to 
the sky of 0. 13° ± 0.02° dGhez et al.ll2008l) . we find a separa- 
tion of p = 588 ± 1 mas and position angle PA = 157.4° ±0.2° 
east of north. If it is bound to KELT-1, it has a projected phys- 
ical separation of ~ 154 ± 8 AU, and a period of ~ 1700 years 
assuming a circular, face-on orbit. 

We used the Galactic model from lDhital et alj (1201 Ot) to as- 
sess the probability that the companion is an unrelated star 
(i.e., a chance alignment). The model uses empirical number 
density distributions to simulate the surface density of stars 
along a given line-of-sight and thus determine probability of 
finding a star within a given angular separation from KELT- 1 . 
We estimate an a priori probability of ~ 0.05% of finding a 
star separated by < 0.59" from KELT- lb. We therefore con- 
clude that the companion is likely to be a bona fide, physi- 
cally associated binary system. With a total proper motion 
of ~ 20 mas/year, it will be possible to definitively determine 
whether the candidate tertiary is physically associated with 
KELT- 1 within one year. 

We note that the companion is unresolved in our follow-up 
primary transit photometry, and thus in principle leads to a 
dilution of the transit signal and a bias in the para meters we 
infer from a fit to the photometry described in § 15.21 However, 
the effect is negligible. As we discuss in the next section, 
we are confident that the primary is being eclipsed. Thus the 
fractional effect on the transit depth is of the same order as the 
fractional contribution of the companion flux to the total flux, 
which is < 1%. 

4. EVIDENCE AGAINST A BLEND SCENARIO 

One of the many challenges of photometric surveys for tran- 
siting planets is the relatively high rate of astrophysical false 
positives, blended eclipsing stellar binary or triple systems 
that can mimic some of the observable signatures of tran- 
siting low-mass companions to single stars. In the case of 
the KELT-North survey, one generic ally expects a higher rate 
of false positives as compared to other wide-field transit sur- 
veys such as HATNet or SuperWASP, because of the poorer 
image quality arising from the comparatively smaller aper- 
ture, larger pixel scale, and wider FOV. For KELT-1 in par- 
ticular, the extreme properties of the companion, relatively 

21 http://www2.keck.hawaii.edu/inst/nirc2/forReDoc/post_observing/dewarp/ 



high vsin/* of the primary, and the fact that the primary is 
somewhat evolved, are all reminiscent of f alse positives that 
have b een reported in previous surveys, e.g., Mandushe v et al.1 
J2005I 

In the case of KELT- lb, however, we have a number of lines 
of evidence that strongly disfavor a blend scenario. 

The most compelling arguments against blend scenarios 
arise from the spectra. First is the lack of strong evidence for 
bisector span variations. The lower panel of Figure |3]shows 
the bisector variations phased to the best-fit period of the com- 
panion as determined from the joint fit to the RV and photom- 
etry data described in §15.21 There is no evidence for bisector 
variations correlated with the orbital phase of the companion. 
The weighted RMS of the bisector spans, excluding the data 
taken on UT 2012-01-07, is - 120 m s" 1 , only - 30% larger 
than would be expected based on the native uncertainties, and 
a factor of ~ 30 times smaller than the RMS of the RV mea- 
surements themselves. Figure [9] shows the bisector spans as 
a function of radial velocity relative to the system barycen- 
ter. There is no strong correlation; the correlation coefficient 
is only -0.17. In contrast, Figure [TOl shows data taken on the 
night of UT 2012-01-07, which covered the primary transit. 
For the subset of these data taken within 0.03 days of the tran- 
sit center (approximately the middle half of the transit), there 
is a clear correlation between the radial velocity and the bi- 
sector variations, with a correlation coefficient of 0.68. This 
is expected since the anomalous radial velocity shift from the 
RM effect is due to a distortion of the rotationally-broadened 
stellar spectral lines as the planet progressively occults the 
light from different parts of the face of the star. Indeed, the 
second piece of evidence that the transit signatures are indeed 
due to a small companion occulting the primary star is the RM 
signal itself (Fig. HJ, which has an amplitude consistent with 
the apparent transit depth and spectroscopically-determined 
vsin/*. 

Third, photometric observations in several different filters 
(n'z) are all consistent with the primary transit having nearly 
the same depth, and are well-modeled by transits of a dark 
companion across a star with the limb darkening consistent 
with its spectroscopically measured T s ff and logg (see Section 
O- 

Fourth, photometric observations at the predicted time 
of superior conjunction reveal no evidence for a secondary 
eclipse at the < 1 mmag level. These first two pieces of 
evidence tend to exclude or strongly disfavor blend sce- 
narios in which the observed transits are due to diluted 
eclipses of a muc h fainter and redder eclipsing binary (e.g., 
lO'Donovan etal]l2006h . 

Finally, our adaptive optics imaging does not reveal any 
sources further than ~ 0.25"from the primary that could be 
both blended with it in seeing-limited images and cause tran- 
sits at the observed depth of ~ 1%. The one source we do 
detect, the putative tertiary, has a flux ratio relative to the pri- 
mary of only ~ 0.5% in the near-IR, and is likely considerably 
fainter in the optical, and thus is too faint to explain the ob- 
served transits. 

We did not perform any detailed modeling to determine 
the viability of sp ecific blend scenarios. We defer here to 
iBakos et"al] d2012l) . who argue that such analyses are gener- 
ally unnecessary in situations in which there are no signif- 
icant bisector variations, the transit ingress/egress durations 
are short compared to the total duration, and the radial veloc- 
ity variations phase with the predicted transit ephemeris. 

We conclude that all of the available data are best explained 
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FIG. 9. — Bisector spans versus the RV relative to the system barycenter, 
excluding observations taken on the night of the primary transit on UT 2012- 
01-07. There is no evidence of a significant correlation between the bisector 
and RV variations, and the RMS of the bisector span variations is ~ 30 times 
smaller than the RMS of the RV measurements. 
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FIG. 10. — bisector spans versus the RV relative to the system barycenter 
for observations taken on the night of the primary transit on UT 2012-01-07. 
The points in red are the subset of those data that were taken within 0.03 days 
of the center of the primary transit, roughly corresponding to the middle half 
of the full transit duration. Note that these data are strongly correlated with 
the RV variations due to the RM effect. 

as due to a Jupiter-sized, brown dwarf companion transiting 
a rapidly-rotating mid-F star, with little or no evidence for 
significant contamination from blended sources. Under this 
assumption, we proceed in the following section to analyze 
these data in order to determine the physical properties of the 
KELT-1 host star and its short-period, low-mass companion. 

5. CHARACTERIZATION OF THE STAR, COMPANION, AND ORBIT 

5.1. Properties of the Host Star 

Table[3]lists various collected properties and measurements 
of the KELT-1 host star. Many these have been culled from 
the literature, and the remainder are derived in this section. In 
summary, KELT- 1 is a mildly evolved, solar-metallicity, mid- 
F star with an age of ~ 1.5-2 Gyr located at a distance of 
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FIG. 1 1 . — The red errorbars indicate measurements of the flux of KELT-1 
host star in various optical and IR passbands. The vertical errorbar indicates 
the photometric uncertainty, whereas the horizontal errorbar indicates the ef- 
fective width of the passban d. The solid curve is th e best-fit theoretical SED 
from the NextGen models of Hauschildt et al. 1 1999), assuming 7" e ff, logg and 
[Fe/H] fixed at the values in Table|4] with A y and d allowed to vary. The blue 
dots are the predicted passband-integrated fluxes of the best-fit theoretical 
SED. 

i~ 260 pc, with kinematics consistent with membership in the 
thin disk. 

We construct an empirical spectral energy distribution 
(SED) of KELT-1 using optical flu xes in the Bj a nd Vj 
passbands from the Tycho-2 catalog (H0g et al. 2000), near- 
infrared (IR) fluxes in the J, H an d Ks passband s from 
the 2MASS Point Source Catalog dSkrutskie et al.l 120061: 
ICutrietal.ll2003l) . and near- and mid-IR fluxes in the four 
WISE passbands (IWright et al.l I20TI ICutri & et al.1 120121) . 
This SE P is shown infill We fi t this SED to NextGen mod- 
els from Hausc hildt et al.l (1 19991) by fixing the values of T e ff, 
logg and [Fe/H] inferred from the global fit to the light curve 
and RV data as described in 35.21 and listed in Table @] and 
then finding the values of the visual extinction Ay and dis- 
tance d that minimizes x 2 . We find Ay = 0.18 ± 0.10 and 
d = 262 ± 14pc, with a X 2 = 10.5 for 6 de grees of freedom, in- 
dicating a reasonable fit (P(> x 2 ) ~ 10%). We also performed 
a fit to the SED without priors, finding T eii = 6500 ± 400K, 
A v = 0.20 ±0.15, logg = 4.25 ± 0.75 and [Fe/H]= -0.5 ± 0.5, 
consistent with the constrained fit. There is no evidence for 
an IR excess. 

We note that the quoted statistical uncertainties on Ay and 
d are likely to be underestimated, because we have not ac- 
counted for the uncertainties in values of T e s, logg and [Fe/H] 
used to derive the model SED. Furthermore, it is likely that al- 
ternate model atmospheres would predict somewhat different 
SEDs and thus values of the extinction and distance. 

In Figure [12] we plot the predicted evolutionary track of 
KELT-1 on a theoretical HR diagram (logg vs. r e ff) , from 
the Yonsei-Yale stellar models (Demarq ue et al.ll2004l) . Here 
again we have u sed the values of M* and [Fe/H] derived from 
the global fit ( 35.21 and Table |4)l. We also show evolution- 
ary tracks for masses corresponding to the ±1 a extrema in 
the estimated uncertainty. In order to estimate the age of 
the KELT-1 system, we compare these tracks to the values 
of r e ff and logg and associated uncertainties as determined 
from the global fit. These intersect the evolutionary track for 
a fairly narrow range of ages near ~ 2 Gyr. The agreement 
between the prediction from the evolutionary track at this age 
and the inferred temperature and surface gravity for KELT- 1 
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FIG. 12. — Theoretical HR di agram based on Yonsei-Yale stellar evolu- 
tion models IDemarque et al. 2004). The solid track is the isochrone for the 
best-fit value s of the mass and metalficity of the host star from the joint fit 
described in g/I] M, = 1.32 ±0.03 M and [Fe/H]= 0.01 ±0.07. The red 
cross shows the best-fit T eft = 6513 ± 50K and logg = 4.229^ from the 
final analysis. The black cross shows the inferred 7" c ff and logg from the 
spectroscopic analysis alone. The blue dots represent the location of the star 
for various ages in Gyr. The host star is slightly evolved with a probable age 
of ~ 2 Gyr, although a similar analysis with a different stellar evolutionary 
model prefers a slightly younger age of ~ 1.5 — 1.75 The two sets of dashed 
tracks represent the tracks for the extreme range of the l-cr uncertainties on 
M, and [Fe/H] from the final analysis (dark shaded), and from the spectro- 
scopic constraints alone (light shaded). 

is remarkably good, but perhaps not entirely surprising. The 
values of T e s, logg, [Fe/H] and M* were a ll de termined in the 
global fit to the light curve and RV data in §15.21 which uses the 
empirical r elations betw e en (T^ ff, logg, [Fe/H]) and (M»,/?„) 
inferred by iTorres et"afl (1201 0b as priors on the fit, in order 
to break the well-known degeneracy between M* and for 
single-lined spectroscopic eclipsing systems. These empiri- 
cal relations are known to reproduce the constraints between 
these parameters imposed by th e physics of stellar evolution 
quite well (see, e.g., Section 8 in Torres et al. 2010|). 

Based on its T e s and J — K and the empirical table of 
spectral type versus c olor a nd T e ff for main-sequence from 
iKenyon & Hartmannl d!995h . we infer the spectral type of 
KELT- 1 to be F5 with an uncertainty of roughly ± 1 spectral 
type. 

We determined the Galactic U,V,W space velocities of 
the KELT-1 system using the proper motion of {pi a ,^s) = 
(-1 1 ±0.7,-9.4 ± 0.7) m as yr" 1 from the NOMAD cata- 
log (IZacharias et alj 120041) . the distance of d = 262 ± 14 pc 
from our SED fit described above, and the barycentric ra- 
dial velocity of the system as determined from the TRES 
observations f ^3.3b of j a bs = -14.2 ±0.2 km/s. We used a 
modification of the IDL routine GAL_UVW, w h ich is itself 
based on the method of iJohnson & Soderbloml d!987l) . We 
adopt the correction for the S un's motion with r e spect to the 
Local Standard of Rest from ICoskunoglu et al.l (1201 ll) . and 
choose a right-handed coordinate system such that positive U 
is toward the Galactic Center. We find (U,V,W) = (19.9 ± 
1.1,-9.6 ± 0.5,-2.6 ± 0.9)km s" 1 , cons istent with member- 
ship in the thin disk (Bensby et al. 2003). We note also that 
the distance of KELT-1 from the Galactic plane is ~ 80 pc. 

Finally, we use the solar evolutionary models of 
Guenther et al.l (fl992|), updated with input physics from 
van Saders & Pinsonneaultl (120121) . to gain some insight into 
the detailed structure of the host star. Fixing the mass 
and metallicity at the values determined from the global fit 



( 35.21 and Tabled, we evolved the model forward until the 
model logg and T e ff approximately matched the values in- 
ferred for KELT-1. We found that ages of ~ 1.5— 1.75 Gyr 
best matched the available constraints, and thus this model 
prefers a so mewhat younger age th an the Yale-Yonsei model 
of (Fig. [12] IDemarque et al.ll2004l) . We therefore decided to 
adopt an age of 1.75 ±0.25 Gyr, consistent with both esti- 
mates. For this range of ages, the models of Guenther et al. 
( 1992) predict a radius of the base of the convective zone of 
R cz = 1.30 ±0.03 Rq, and a very small mass for the convec- 
tive zone of M cz = [2.8 ±0.14] x 10" 5 M , as expected given 
the effective temperature of r e ff ~ 6500K. In addition, the 
moment of inertia for the star and convective zone are C* = 
[1. 15 ±0.04] x 10 54 gcm 2 andC fz = [3.2±0.6] x 10 51 gcm 2 , 
respectively. We can also write the m oment of inertia of th e 
star as C* = a*M*R 2 * with a» = 0.0422 dGuenther et al.lll992l) . 
We will use these to estimate the angular momenta of the star, 
companion and orbit in 36.21 

5.2. System Properties Derived from a Joint Fit 

It is well known that a joint fit to high-quality RVs and 
transit photometry of a transiting planet system allows one 
to determine the mass and radius of the star and planet, 
as well as the semi-major axis of the orbit, in principle to 
very high precision, up to a perfect one-parameter degener- 
acy dSeager & Mailen-OrnelasH2003l) . This degeneracy arises 
because the duration, depth, and shape of the primary tran- 
sit, when combined with the eccentricity, longitude of pe- 
riastron, and period of the planet from RV data, allow one 
to precisely estimate the density of the primary star but 
not M.J, or R* separately. Breaking this —R* degeneracy 
generally requires imposing some e xternal constraint, such 
as a theoretical mass-radius rela tion dCodv & Sas selovl 120021 : 
Seager & Mallen-Ornelas 200 3]), or con straints from theoret- 
ical isochrones (e.g. jBakos et al.ll2012l) . In principle, a mea- 
surement of logg from a high-resolution spectrum can be used 
to break the degeneracy, but in practice these measurements 
are generally not competitive with the constraint on and 
often have systematic uncertainties that are comparable to the 
statistical uncertainties. 

We fitted the RV and transit data us ing a modified version of 
the ID L fitting package EXOFAST dEastman, Gaudi & Agoll 
120121) . The approach of EXOFAST to breaking the M* - 
R* degeneracy is similar to the method described in, 
e.g.. I Anderson et all d2012l) . but with significant differences. 
We will review it briefly he re, but point the reader to 
Eastman, Gaudi & Agol (2012) for more details. We fitted the 
RV and tran sit data simultaneously with a standard Keplerian 
and transit (Mandel & Agol 2002) model using a modified 
MCMC method (described in more detail below). In addition 
to the standard fitting parameters, we also included r e ft,logg, 
[Fe/H] as proposal parameters. We then included priors on 
the measured values of r e ff,logg, [Fe/H] as determined from 
analysis of the TRES spectra and given in 33.31 In addition, 
we included priors on M» and which are based on the 
empirical relati ons between (T^ffjogg, [Fe/H]) and (M,,/?*) 
determined by ITorres et al.l (|2010). These priors essentially 
break the M*-R* degeneracy, as they provide similar con- 
straints as isochrones, i.e., they encode the mapping between 
the M*, [Fe/H] and age of a star to its T e g and logg as dictated 
by stellar physics. 

We fitted the 6 primary transits, Doppler RV, stellar param- 
eters, and RM effect simultaneously using EXOFAST, which 
employs a Differential Evolution Markov Chain Monte Carlo 
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(DE-MC) method dterBraakll2006l). We conve rted all times 
to the BJDtdb standard (lEastman et al.l 1201 Obi) , and then at 
each step in the Markov Chain, we converted them to the tar- 
get's barycentric coordinate system (ignoring relativistic ef- 
fects). Note the final times were converted back to BJDtdb 
for ease of use. This transformation accurately and transpar- 
ently handles the light travel time difference between the RVs 
and transits. 

First, we fitted the Doppler RV data independently to a sim- 
ple Keplerian model, ignoring the RM data taken on UT 2012- 
01-07. At this stage, we did not include any priors on the 
stellar parameters, as they do not affect the RV-only fit. We 
scaled the uncertainties such that the probability that the \ 2 
was larger than the value we achieved, P (> % 2 J, was 0.5, to 
ensure the resulting parameter uncertainties were roughly ac- 
curate. For a uniform prior in eccentricity, we found the orbit 
is consistent with circular, with a 3a upper limit of e < 0.04. 
Nevertheless, in order to provide conservative estimates for 
the fit parameters, we allowed for a non-zero eccentricity in 
our final fit. However, to test the effect of this assumption, we 
repeated the final fit forcing e = 0. We also investigated the 
possibility of a slope in the RVs, but found it to be consistent 
with zero, so we did not include this in the final fit. 

Next, we fitted each of the 4 transits individually, includ- 
ing a zero point, Fqj and airmass detrending variable, Co,; for 
each of the i transits. The airmass detrending coefficient was 
significant (> lcr) for all but 1 transit, so for consistency, we 
included it for all. Aft er finding the best fit with AMOEBA 
dNelder & MeadlH%5h . we scaled the errors for each transit 
such that P (> x 2 ) = 0.5. At this stage, we included the priors 
on the stellar parameters as described above. 

Next, we performed a combined fit to all the data, includ- 
ing a prior on the projected stellar rotation velocity ( vsin/* = 
55.2 ±2 km s -1 ) from the spectra 22 , and a prior on the period 
from the KELT-North discovery light curve (P = 1.217513 ± 
0.000015 days). Because it is usually systematics in the RV 
data that vary over long time scales (due to a combination of 
instrumental drift and stellar jitter) that ultimately set the er- 
ror floor to the RVs, we expect the uncertainties of densely 
packed observations to be smaller than the RMS of all ob- 
servations, but with a systematic offset relative to the rest of 
the orbit. Therefore, we fitted a separate zero point during 
the RM run, and also scaled the errors on the RVs during 
transit to force P (> x 2 ) = 0.5 for those subsets of points. 
Pi> x 2 ) depends on the number of degrees of freedom, but 
it is not obvious how many degrees of freedom there are in 
the RM run - technically, the entire orbit and the transit af- 
fect the x 2 °f m e RM (13 parameters), but the freedom of the 
RM measurements to influence most of those parameters is 
very limited, when fit simultaneously with the transits. In- 
deed, even vsin/* is constrained more by the spectroscopic 
prior than the RM in this case, which means there are only 
two parameters (the projected spin-orbit alignment, A, and the 
zero point, 7) that are truly free. To be conservative, we sub- 
tracted another degree of freedom to encompass all the other 
parameters on which the RM data has a slight influence, be- 
fore scaling the errors. 

The RM data were modeled using the lOhta et all (120051) 
analytic approximation with linear limb darkening. At each 

22 The prior on vsin/<. improves t he determination of the spin-orbit align- 
ment angle A (Gaudi & Winn 2006). We also performed a fit without this 
prior, finding results that were roughly consistent with, although less precise 
than, those with the prior. 



step in the Markov Chain, we interpo l ated t he linear limb 
darkening tables of Claret & Bloe menl (1201 lb based on the 
chain's value for logg, J e ff7 and [Fe/H] to derive the linear 
limb-darkening coefficient, u. We assumed the V band pa- 
rameters to approximate the bandpass of TRES, though we 
repeated the exercise in the B-band with no appreciable dif- 
ference in any of the final parameters. Note that we do not fit 
for the limb-darkening parameters, as the data are not suffi- 
cient to constrain them directly.. The uncertainties in all the 
limb-darkening parameters provided in Table arise solely 
from the scatter in logg, T e ff, [Fe/H]. We assume no error in 
the interpolation of the limb-darkening tables. 

In order to search for Transit Timing Variations (TTVs), 
during the combined fit, we fitted a new time of transit, 7c., 
for each of the i transits. Therefore, the constraint on Tc and 
P (quoted in Tables|5]and|4l respectively) come from the prior 
imposed from the KELT-North light curve and the RV data, 
not the follow up light curves. Using these times to constrain 
the period during the fit would artificially reduce any observed 
TTV signal. A separate constraint on Tc and P follows from 
fitting a linear ephemeris to the transit times, as discussed in 
35.31 It is the result from this fit that we quote as our final 
adopted ephemeris. 

The results from this global fit are summarized in Tables |4] 
and [5] We also show the results for the physical parameters 
assuming e = in Table HI the differences between the fixed 
and free eccentricity fits are always smaller than their uncer- 
tainties, and generally negligible for most of these parame- 
ters. The values of 7" e ff, logg, and [Fe/H] we infer from the 
global fit are in agreement with the values measured directly 
from the TRES spectra to within the uncertainties. Since the 
spectroscopic values were used as priors in the global fit, this 
generally indicates that there is no tension between the value 
of p* inferred from the light curve and RV data, and the spec- 
troscopic values. The median value and uncertainty for T s g 
is nearly unaffected by the global fit. While the median value 
for [Fe/H] has changed slightly, the uncertainty is very similar 
to that from the input prior. On the other hand, the uncertainty 
in logg from the global fit is a factor of > 5 smaller than the 
uncertainty from the spectroscopic measurement. This is not 
surprising, since the constraint on p* from the RV and light 
curve d ata provides a stro ng constraint on log g via the rela- 
tions of IsSsetaDdZQlS). 

Following papers by the HATNet collaboration (e.g., 
IBakos et al.ll201 U IHartman et al.l 1201 lb . we also present in 
Table |4] our estimates of the median values and uncertainties 
for a number of derived physical parameters of the system 
that may be of interest, including the planet equilibrium tem- 
perature assuming zero albedo and perfect redistribution T eq , 
the average amount of stellar flux at the orbit of the compan- 
ion (F), and the Safrono v number 6 = (a/R p )(Mp/M*) (e.g., 
Hansen & Barman 2007). In addition, in Table [5] we quote 
our estimates of various fit parameters and intermediate de- 
rived quantities for the Keplerian RV fit, the primary transits, 
and the secondary eclipse. 

We note that the final uncertainties we derive for M„, 
logg and p* are relatively small, ~ 2% -5%. These uncer- 
tainties are similar to those found for other transiting planet 
systems using metho ds similar to the one used here (e.g., 
lAnderson et al.l 120121) . Specifically, these methods derive 
physical parameters from a fit to the light curve and RV data, 
which simultaneously imposes an empirical constraint be- 
twee n the M t , R* and r e ff , log and [Fe/H] ultimately derived 
from iTorres et al.l (120101) . This constraint helps to break the 
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M* —R * degeneracy pointed out by Seager & Mal len-Ornelasl 
(2003) and discussed above, and ultimately dictates our fi- 
nal uncertainty on and (and thus Mp and Rp). In 
particular, our spectroscopic measurement of log g provides 
a much we aker constraint. Giv en that our results rely so heav- 
ily on the Torr es et all (|2010) relations, it is worthwhile to 
ask to what extent our parameters and uncertainties might 
be affected should these relations be systematically in error. 
First, as already noted, these empirical relations are known 
to agr ee well with stellar isochrones in general dTorres et al.l 
2010), and for KELT-1 in particular (Fig. [T2l. Second, anal- 
yses using stellar isochrones rather than empiri cal relations 
produ ce similar uncertainties on M* and R* (e.g.. lBakos et alJ 
1201 lb . suggesting that the small uncertainties w e derive are 
not a by-product of our methodology. Finally, Southworth 
(2009) demonstrated that the results of the analysis of 14 tran- 
siting systems with several different sets of isochrones gener- 
ally agree to within a few percent. We therefore conclude that 
our results are likely accurate, with systematic errors at the 
level of our statistical uncertainties (a few percent). 

5.3. System Ephemeris and Transit Timing Variations 

Table [6] lists the measured transit times for each of the six 
modeled transits, and Figure Q~3] shows the residuals of these 
times from a fit to a linear ephemeris. The best fit has 



T C (B JDtdb) = 2455909.292797 ± 0.00024 
P= 1.2175007 ±0.000018, 



(6) 



which is consistent with the ephemeris derived from the 
KELT-North light curve alone. The \ 2 = 44.9 for the linear 
fit with 4 degrees of freedom is formally quite poor. This 
is mostly driven by one nominally significant (5c) outlier, 
specifically for the transit observed on UT 2011-12-16 from 
FLWO. We note that the faint companion to KELT-1, if indeed 
bound, is too distant to explain such large TTVs. 

We have taken great care to ensure the accuracy of our 
clocks and our conversion, and the fact that the residuals 
from different observatories roughly follow the same trend 
in time suggests that a catastrophic error in the observatory 
clock cannot be blamed. Since we fit the trend with airmass 
simultaneously with the transit, the potentially-large covari- 
ance between it and the transit time should be accurately re- 
flecte d in the quoted uncertainties ( Eastman. G audi & Agoll 
2012). Nevertheless, ourMCMC analysis does not adequately 
take into account the effect of systematic uncertainties, and 
in particular we do not account for cor related uncertainties 
dPont et alJl2006ah ICarter & Winnll2009l) . which could skew 
the transit time of a given event substantially. And, given 
the results from Kepl er which suggest th e rarity of such tran- 
sit timing variations (Stef fen et alJl20!2h . we are reluctant to 
over-interpret this result. Nevertheless, this is an interesting 
target for future follow-up. 

5.4. Secondary Eclipse Limits 

We observed the predicted secondary eclipses of KELT- lb 
assuming e = on UT 201 1-12-02, 201 1-12-30, and 2012-01- 
04. In none of the three were we able to detect a secondary 
eclipse. The observations on 2011-12-02 and on 2012-01-04 
were taken from the ULMO Observatory in i. Both nights 
we were able to observe through the predicted ingress and 
egress of the potential secondary. The two i band light curves 
have a combined 236 data points, and show an RMS scatter of 
1.56 mmag. The observations on 201 1-12-30 were taken with 
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FIG. 13. — The residuals of the transit times from the best-fit (linear) 
ephemeris. The transit times are given in Table|6] 
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FIG. 14. — Values of the heat redistribution parameter /' and Bond albedo 
Ag that are excluded at a given confidence level based on the data taken dur- 
ing the secondary eclipse shown in Figure [6] The /' parameter describes the 
efficiency of heat redistribution, and is 1/4 in the case of uniform redistribu- 
tion, and 2/3 in the case of no redistribution. In between these two extremes 
/' is not easily related to the amount of heat redistribution. The orange sec- 
tion corresponds to eclipse depths detectable at less than the 68% confidence 
level, the light orange for 68-90% confidence, yellow for 90-95% confidence, 
and the white for > 95% confidence. 

the FTN telescope in Pan-STARRS-Z. In this case we were 
only able to begin observations half way through the predicted 
secondary eclipse. This Z-band light curve has 72 points, and 
an RMS scatter of 0.75 mmag. 

We used the system parameters derived from the joint fit 
( i35.2| > to fit our three observations. Since we did not detect 
a secondary eclipse, we used these fits to explore the combi- 
nation of heat redistribution efficiency and Bond albedo As 
that would give rise to a secondary eclipse depth that is in- 
consistent with our data. To do so, we calculated the sec- 
ondary eclipse depths we would expect for a range of redis- 
tribution efficiencies and albedos, and then fit a secondary 
eclipse model with the predicted depth to all three of our ob- 
servations simultaneously. 

In calculating the expected secondary eclipse depths, we 
made the assumption that both the star and the planet were 
blackbodies. We also assumed that the planet was a grey Lam- 
bert sphere, so the geometric albedo A g = (2/3)A#, and the 
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spheri cal albedo is co nstant as a function of wavelength. Fol- 
lowing Seager (2010), we parametrized the heat redistribution 
efficiency as /', which is 1/4 in the case of uniform redistri- 
bution, and 2/3 in the case of no redistribution. We note that 
in between these two extremes /' is not easily related to the 
amount of heat redistribution, i.e., /' = 0.45 does not imply 
half of the incident stellar energy is redistributed around the 
planet. 

To test these expected secondary eclipse depths against our 
observations, we fit simple trapezoidal eclipse curves with the 
expected depths to all three datasets simultaneously. Under 
our assumptions, the depth, timing, shape, and duration of the 
secondary eclipse are all determined by the parameters de- 
rived from the global fit and our specified values of Ag and /'. 
We then fit this model to our data, allowing for a normaliza- 
tion and a linear trend in the flux with time. We used the A% 2 
between the best fit eclipse model and the best constant fit, 
which itself was allowed a free slope and offset, to evaluate 
the detectability of each of the secondary eclipse depths. We 
used the \ 2 distribution to transform these A\ 2 values into de- 
tection probabilities. Figure [6] shows an example light curve 
against a median binned version of our data. This particular 
curve is the secondary eclipse we would expect if KELT- lb 
had Ag = 0. 1 , and instantaneously re-radiated its incident stel- 
lar flux, i.e., /' =2/3. We would have detected this event with 
more than 95% confidence. 

Figure [14] shows the results of our exploration of the heat 
redistribution versus Bond albedo parameter space. The or- 
ange section corresponds to eclipse depths detectable at less 
than the 68% confidence, the light orange is for depths de- 
tectable with 68-90% confidence, light yellow is for 90-95%, 
and the white contains depths that would have been detected 
at greater than 95% confidence. The particular shapes of the 
contours on this plot come from the competing effects of re- 
flection and blackbody emission from KELT- lb on the depth 
of the secondary eclipse. Along the very bottom of the figure 
the Bond albedo is zero, and thus there is only thermal emis- 
sion. We see the strong change in eclipse depth as amount 
of heat redistribution decreases, thus causing the temperature 
and eclipse depth for KELT- lb to increase. Along the top 
of the figure, where the Bond albedo is 0.75, the reflected 
starlight dominates the blackbody emission such that chang- 
ing the redistribution efficiency has little effect on the eclipse 
depth. 

Slightly more than half of the allowed parameter space in 
Figure [14] would have caused secondary eclipses detectable 
in our data at greater than 90% confidence, while almost all 
would have been detected at more than 68% confidence. Since 
we did not see a secondary eclipse in our observations, we 
conclude that KELT- lb either it has a non-zero albedo, or it 
must redistribute some heat from the day side, or both. For- 
mally, the scenario that is most consistent with our data is that 
KELT- lb has both a low Bond albedo and is very efficient at 
redistributing heat away from its day side, however we are 
reluctant to draw any strong conclusions based on these data. 

6. DISCUSSION 

From our global fit to the light curves and RVs, we find 
that KELT- lb is a low-mass companion with a measured mass 
M P 

on a circular orbit with a semimajor axis of a = 0.02466 ± 
0.00016AU. The host KELT-1 is a mildly evolved mid- 
F star with a mass M* = 1.324 ±0.026 M Q , radius R« = 
1.462^;^ R Q , effective temperature T e{f = 6518 ±50 K, and 
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FIG. 15. — (Top panel) Mass versus period for the known transiting 
companions to main-sequence stars with companion masses in the range 
1 — 100 Mj. An estimate of the deuterium burning limit (Spiegel et al. 201 1) 
is shown as the horizontal dotted line, while the hydrogen burning limit is 
shown as the horizontal dashed line. The vertical line shows the divisio n 
between hot and cool stars of T e ff = 6250K suggested by Winn et al. 1 2010). 
Brown dwarfs are shown as triangles, exoplanets as squares, and low-mass 
stars as asterisks. KELT-lb is shown as the large star. It is the shortest period 
transiting brown dwarf currently known. (Bottom panel) Mass versus host 
star effective temperature T e ff for the sample of transiting companions shown 
in the top panel. As suggested by Bouchy et al. (201 lb), there is some evi- 
dence that massive (Mp > 5 Mj up companions are preferentially found around 
hot (Ttff > 6000K) stars, and KELT-lb follows th is possible trend. Note that 
we exclude the BD companion to NLT T 41135 f lrwin et al. 2010), and the 
double BD transit system 2M0535— 05 I Stassun et al. 2007) in this and sub- 
sequent plots. 

a likely age of ~ 1.5-2Gyr. Because of its small semima- 
jor axis and hot host, KELT-lb receives a large stellar inso- 
27.23!^ M Ju P and radius' R P = 1.110ig;g|| R Jup . It is lation flux of {F} = 7.81^ x 10 9 erg s" 1 cm" 2 , implying a 

high equilibrium temperature of T eq = 2422^ K assuming 
zero albedo and perfect redistribution. Both the surface grav- 
ity and density of KELT-lb are substantially higher than that 
of its host star, and higher than we would expect for a stellar 
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object. We find that the orbit normal of KELT- lb is well- 
aligned with the projected rotation axis of its host star, with a 
projected alignment angle of A = 2 ± 16 degrees. 

Even among the large and diverse menagerie of known 
transiting exoplanets and low-mass companions, KELT- lb is 
unique. First, it is one of only 7 unambiguous objects with 
mass the range ~ 13-80 Mj up that are known to transit stars. 
Among these, it has the shortest period, and orbits the bright- 
est host star (V = 10.7). In addition, there is potentially a stel- 
lar M dwarf companion to the primary. For all these reasons, 
KELT- lb is likely to be a very interesting object for further 
study, and we expect it will provide a benchmark system to 
test theories of the emplacement and evolution of short pe- 
riod companions, as well the physics of tidal dissipation and 
irradiated atmospheres of substellar objects. We will discuss 
some of these ideas briefly. 

6.1. Brown Dwarf or Supermassive Planet? KELT-lb and 
the Brown Dwarf Desert 

Is KELT-lb a brown dwarf (BD) or a is it a suppermas- 
sive planet? By IAU convention, brown dwarfs (BDs) are 
defined to hav e masses between t he deuterium burning limit 
of ~ 13 Mj up (Spie gel et al.ll201 1[) and the hy drogen burning 
limit of ~ 80 Mt, t (e.g.. iBurrows et al.lll997l) . Less massive 
objects are defined to be planets, whereas more massive ob- 
jects are stars. By this definition, KELT-lb is a low-mass 
BD. However, it is interesting to ask whether or not KELT- 
lb could have plausibly formed in a protoplanetary disk, and 
therefore might be more appropriat e considered a "supermas- 
sive planet" ( [Schneider et al.ll20Tlh . More generally, it is in- 
teresting to consider what KELT-lb and systems like it may 
tell us about the formation mechanisms of close companions 
with masses in the range of 10- 100M Jup . 

One of the first results to emerge from precision Doppler 
searches for exoplanets is the existence of a BD desert, an 
apparent paucity brown dwarf companions to FGK stars with 
periods less than a few years, relative to the f requency of stel- 
lar co mpanions in the same range of periods (Ma rcv & Butlerl 
2000). Subsequent studies uncovered planetary compan- 
ions to such stars in this range of periods in abundance 
dCumm ing et al. 2008), indicating that the BD desert is a lo- 
cal nadir in the mass function of companions to FGK stars. 
The simplest interpretation is that this is the gap between the 
largest objects that can form in a protoplanetary disk, and the 
smallest objects that can directly collapse or fragment to form 
a self-gravitating object in the vicinity of a more massive pro- 
tostar. Therefore, the location of KELT-lb with respect to the 
minimum of the brown dwarf mass function might plausibly 
provide a clue to its origin. 

6.1.1. Comparison Sample of Transiting Exoplanets, Brown 
Dwarfs, and Low-mass Stellar Companions 

In order to place the parameters of KELT-lb in context, we 
construct a sample of transiting exoplanets, BDs, and low- 
mass stellar companions to main sequence stars. We focus 
only on transiting objects, which have the advantage that 
both the mass and radius of the companions are precisely 
known 23 . We collect the transiting exoplanet systems from the 

23 In contrast, for companions detected only via RVs, only the minimum 
mass is known. Of course, one can make an estimate of the posterior proba- 
bility distribution of the true mass given a measure d minimum m ass by adopt- 
ing a prior for the distribution of inclinations (e.g., Lee et al. 201 1). However, 
this procedure can be particularly misleading in the case of BDs: if BDs are 



Exoplanet Data Explorer (exoplanets.org, Wright et al.l201 ll) . 
discarding systems for which the planet mass is not mea- 
sured. We supplement this list with known transiting brown 
dwarfs dDeleuil et al.l20"08tlJohnson et al.l201 UlBouchv et all 
1201 1 alibi : lAnderson et al.ll201 lh. We do not include the sys- 
tem discovered by (llrwin et al .11201 Oh . because a radius mea- 
surement for the brown dwarf was not possib le. We also do 
not include 2M0535-05 (IStassun et al.ll2007l) . because it is a 
young, double BD system. We add sever al transiting low- 
mass s t ars near the hydro gen burning limit ( Pont et al.ll2005[ 
2006b; iBeattv et alj 120071) We adopt the mass of XO-3b 
from the discovery paper dJohns-Krull et al.ll2.Q08h . which is 
Mp = 13.1 zfc0.4 Mjnp, wh ich straddles the deuterium burning 
limit dSpiegel et alj|201lh . However, later estim ates revised 
the m ass significantly lower to Mp =11.8 ±0.6 (Wi nn et ail 
2008). We will therefore categorize XO-3b as an exoplanet. 

The disadvantage of using samples culled from transit sur- 
veys is that the sample size is much smaller, and transit 
surveys ha ve large and generally unquantified selection bi- 
ases (e.g., Gaudi et al. 2005; Fressin et al. 2009), particularly 
ground-based transit surveys. We emphasize that such biases 
are almost certainly present in the sample we construct. We 
have therefore made no effort to be complete. The compar- 
isons and suggestions we make based on this sample should 
not be considered definitive, but rather suggestive. 

Figure Q3] places KELT-lb among the demographics of 
known transiting companions to main sequence stars, focus- 
ing on massive exoplanets, BDs, and low-mass stars with 
short periods of < 30 days. KELT-lb has the tenth short- 
est period of any transiting exoplanet or BD known. It has 
the sixth shortest period among giant {Mp > 0.1 Mj up ) plan- 
ets, with only WASP- 19b, WASP-43b, WASP- 18b, WASP- 
12b, OGLE-TR-56b, and HAT-P-23b having shorter periods. 
KELT-lb is more massive by a factor ~ 3 tha n the most mas- 
sive of these, WASP- 18b dHellier et al.ll2009t) . KELT-lb has a 
significantly shorter period than any of the previously known 
transiting brown dwarfs, by a factor of > 3. KELT-lb there- 
fore appears to be located in a heretofore relatively unpopu- 
lated region of the Mp-P parameter space for transiting com- 
panions. 

Although the KELT-1 system is relatively unique, it is 
worth asking if there are any other known systems that 
bear some resemblance to it. The Mp ~ 18Mj up , P ~ 1.3 
day RV-discoyered companion to the M dwarf HD 41004B 
dZucker et al.ll2003l) has similar minimum mass and orbit as 
KELT-lb, however the host star is obviously quite different. 
Considering the host star prop erties as well, perh aps the clos- 
est analogs are WASP- 18b dHellier et al.l [Iffl) , WASP-33b 
dCollier Cameron et all 120101) and KOI- 13b dMazeh et alj 
l2012t iMislis & Hodgkinll2012l) . All three of these systems 
consist of relatively massive (M p > 3 Mj up ) planets in short 
(< 2 day) orbits around hot (T e ff > 6500 K) stars. 

The mass of KELT-lb (~ 27 Mj up ) is close to the most 
arid part of the BD d esert, estimated to be at a m ass of 
31^8 Mj up according to Grether & Line weaver! d2006 ) . Thus, 
under the assumption that the BD desert reflects the diffi- 
culty of forming objects with this mass close to the parent 
star under any formation scenario, KELT-lb may provide an 
interesting case to test these various models. For disk scenar- 

indeed very rare, then objects with minimum mass in the BD desert are more 
likely to be stellar companions seen at low inclination. Anecdotally, in those 
cases where constraints on the inclinations can be made, companions with 
minimum mass near the middle of the brown dwarf desert often do turn out 
to be stars (e.g., Sa hlmann et alJ201 1 ; Fleming et al.l2012h . 
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ios, gravitational instability can likely form su ch massive ob- 
jects, but likely only at much larger distances dRafikovll200"5l : 
iDodson-Robinson et al.ll2009t iKratteret alj|2010h . The maxi- 
mum mass possible from core ac cretion is poorly unde rstood, 
but may be as large as ~ 40 Mj up dMordasini et al.ll2009l) . The 
possibility of significant migration of KELT- lb from its birth 
location to its present position must also be considered, par- 
ticularly given the existence of a possible stellar companion 
to KELT-1 ( 33.51 ). This possibility complicates the interpreta- 
tion of the formation of KELT- lb significantly. For example, 
it has been suggested that bro wn dwarf companions are mor e 
common at larger separations (Metchev & Hillenbrand 2009); 
thus KELT- lb may have formed by collapse or fragmentation 
at a large separation, and subsequently mi grated to its current 
positi on via the Kozai-Lidov mechanism ( lKozai|[T96l : iLidovl 
fl962h . 

One clue to the origin of KELT- lb and the BD desert may 
be found by studying the frequency of close BD companions 
to stars as a function of the stellar mass or temperature. Figure 
rT~5l shows the mass of known transiting short period compan- 
ions as a function of the effective t e mpera ture of the host stars. 
As pointed out by iBouchv et alj d201 lbl) . companions with 
Mp > 5 Mj up appear to be preferentially found around hot stars 
with 7/eff > 6000 K, and KELT- lb follows this trend. Although 
these hot stars are somewhat more massive, the most dra- 
matic difference between stars hotter and coole r than 6000 K 
is the d epth of their convection zones. This led lBouchy et al.1 
d201 lbl) to suggest that tides may play an important role in 
shaping the frequency and distribution of massive exoplanet 
and brown dwarf companions to old stars. So me evidence 
for this has been reported by IWinn et ail (1201 0b . who argue 
that hot (T e ff > 6250K) stars with close companions prefer- 
entially have high obliquities, suggesting that if the emplace- 
ment mechanisms are similar for all stars, tidal forces must 
later serve to preferentially bring cool host stars into align- 
ment. Figure [16] shows the distribution of spin-orbit align- 
ments for transiting planets versus the host star effective tem- 
perature. KELT- lb falls in the group of hot stars with small 
obliquities. Interestingly the other massive > 5 Mj up planets 
are also located in this group. 

We discuss the possible formation and evolutionary history 
of KELT- lb, and the likely role of tides in this history, in more 
detail below. We remain agnostic about the classification of 
KELT- lb as a brown dwarf or supermassive planet. 

6.2. Tides, Synchronization, and Kozai Emplacement 

Given the relatively large mass and short orbital period of 
KELT- lb, it seems probable that tides have strongly influ- 
enced the past evolution of the system, and may continue to 
be affecting its evolution. The literature o n the influence of 
tides on exoplanet systems is vast (see, e . gjRasio et al| 199a ; 
Qgilvie & Lin1l2q04ilJackson et all 2008; Le conte et alJl2010b 
Matsumura et al.l2010tlHansenl2010l for but a few examples), 
and there appears to be little consensus on the correct treat- 
ment of even the most basic physics of tidal dissipation, with 
the primary uncertainties related to where, how, and on what 
time scale the tidal energy is dissipated. 

While we are interested in evaluating the importance of 
tides on the evolution of the orbit of KELT- lb and the spin of 
KELT-1, delving into the rich but difficult subject of tides is 
beyond the scope of this paper. We therefore take a somewhat 
heuristic approach. Specifically, we construct a few dimen- 
sionless quantities that likely incorporate the primary physi- 
cal properties of binary systems that determine the scale of 
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FIG. 16. — The projected spin-orbit alignment angle A for transiting plan- 
ets as measured by the RM effect, versus the effective temperature of the 
host star, following Winn et al. (2010). The grey points show exoplanets with 
mass Mp < 5 Mj up , whereas the black points show those with Mp > 5Mj up . 
KELT- lb, shown with a star, is the first transiting brown dwarf with a RM 
measurement. Its orbit normal is consistent with being aligned with the pro- 
jected host star spin axis. The dotted vertical line shows the suggested divid- 
ing line between hot and cool stars by Winn et al. (2010). 

tidal evolution, but do not depend on the uncertain physics of 
energy dissipation. Specifically, we define, 



Ta: 



Mp 
~Mp 



and 



(7) 



(8) 



For some classes of theories of tidal dissipation and un- 
der some assumptions, T a is proportional to the e-folding 
timescale for decay of the orbit, and Tui, is proportional to 
the timescale for synchronization of the spin of the star with 
the companion orbital period. It is worthwhile to note that 
for transiting planet systems the combinations of parameters 
Mp/Mt, and a/R* are generally much better determined than 
the individual parameters. In particular, the ratio of the mass 
of the planet to that of the star is closely related to the RV 
semi-amplitude K, whereas a/R* is closely related to th e ratio 
of the duration of the transit to the period dWinnll20T()|) , Fig- 
ure[l7]shows T a and 7w, as a function of orbital period for the 
sample of transiting exoplanets, brown dwarfs, and low-mass 
stars discussed previously. KELT- lb has shorter timescales 
than nearly the entire sample of systems, with the exception 
of a few of the low-mass stars. We therefore expect tidal ef- 
fects to be quite important in this system. 

As a specific example, under the constant time lag model 
dHutl 11981b iMatsumura eTaT] |201Q|) . and assuming dissipa- 
tion in the star, zero eccentricity, zero stellar obliquity, and 
a slowly rotating star, the characteristic time scale for orbital 
decay due to tides is Td ecay = a/\d\ = (127r) _1 Q'^T a P ', where 
Q'^ is related to the dimensionless tidal quality factor. For 
KELT- lb, T a ~ 3 x 10 4 , and so r deca , ~ 0.3 Gyr for = 10 8 , 
clearly much shorter than the age of the system. Similarly, 
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FIG. 17. — Dimensionless combinations of physical parameters that quan- 
tify the relative time scale for orbital tidal decay (top panel) and stellar spin- 
orbit synchronization (bottom panel) for differe nt bin ary systems, as a func- 
tion of the orbital period of the system. See 36.21 foi an explanation and 
assumptions. Brown dwarfs are shown as triangles, exoplanets as squares, 
and low-mass stars as asterisks. KELT-lb is shown as the large star. Among 
known transiting exoplanets and brown dwarfs, it has the shortest character- 
istic time scale for orbital decay and synchronization. 

the time scale for spinning u p the star by the companion is 
Tsynch = w*/|ci*| oc Q'^Tuj^P (fMatsumur a et al.ll2010h . and so 
is also expected to be short compared to the age of the system. 

Given the expected short synchronization time scale and 
the fact that the expected time scale for tidal decay is shorter 
than the age of the system, it is interesting to ask whether 
or not the system has achieved full synchronization, thus en- 
suring the stability of KELT-lb. The measured projected ro- 
tation velocity of the star is vsini* = 56 ± 2 km s" 1 , which 
given the inferred stellar radius corresponds to a rotation pe- 
riod of P» = 2izR* sin/*/vsin/* = [1.322±0.053]sin/„ days, 
which differs from the orbital period of KELT-lb by ~ 2a 
for = 90°. This is suggestive that the system is indeed 
synchronized. The small discrepancy could either be due to 



a slightly underestimated uncertainty on vsin/*, or the host 
could be moderately inclined by /» ~ [67 ±7]°. However, 
one might expect the obliquity of the star to be realigned on 
roughly the same time s cale as the synchronization of its spin 
(Matsumur a et al.ll2010b . The stellar inclination can also be 
constrained by the precise shape of the transit light curve: 
lower inclinations imply higher rotation velo cities, and thu s 
increased oblateness and gravity brightening dBarnesI 2009). 
Ultimately, the inclination is limited to > 10° in order to 
avoid break up. 

We can also ask, given the known system parameters, if the 
system is theoretically expected to be able to achie ve a sta- 
ble synchronous state. A system is "Darwin stable" (iDarwinl 
1 1 897t iHutlfl 980h if its total angular momentum, 



Ltot —L or b+ L u i * + L M p 
is more than the critical angular momentum of 
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where L orn is the orbital angular momentum, L w>> is the 
spin angular momentum of the star, L w p is the spin an- 
gular momentum of the planet, and C* = a*M*Rl and 
Cp = apMpRp are t he moments of inertia of the star and 
planet, respectively (Matsumura et al. 2010). Since Cp/C* ~ 
(Mp/M«)(R P /R*) 2 ~ 10" 3 , the contribution from the planet 
spin to the total angular momentum is negligible. We find 
Ltot /L cr it = 1 .029 ± 0.014, marginally above the critical value 
for stability. In addition, we find (L u + L w p) / 'L or b = 0.154 ± 
0.006, which is smaller than th e maxim um value of 1/3 re- 
quired for a stable equilibrium (lHutlll980l) . Curiously, if we 
assume the star is already tidally synchronized, we instead in- 
fer (L Wi *+L Wi p)/L or j = 0.167 ± 0.004, i.e., remarkably close 
to exactly one-half the critical value of 1/3. 

Two additional pieces of information potentially provide 
clues to the evolutionary history of this system: the detection 
of a possible tertiary ( jj3.5l FigJS), and the measurement of 
the RM effect (Fig.|4]i, demonstrating that KELT-1 has small 
projected obliquity. If the nearby companion to KELT-1 is 
indeed bound, it could provide a way of emplacing KELT- 
lb in a small orb it via the Kozai-Lidov mechanism dKozail 
[T96llLidovli962h . If KELT-lb were originally formed much 
further from its host star, and on an orbit that was signifi- 
cantly misaligned with that of the putative tertiary, then its 
orbit might subsequently be driven to high eccentricity via 
secular perturbations from the tertiary dHolman et all 11997b 
iLithwick & Naozll201 U iKatz et alj|201 lb . If it reached suffi- 
ciently high eccentricity such that tidal effects became impor- 
tant at periastron, the orbit wo uld be subsequently circular- 
ized at a relatively short period (jFabrvckv & Tremainell2007[ 
IWu et alJl2007l: ISocrates et al.ll2012l) . Nominally, one might 
expect the orbit of KELT- lb to b e then left with a relatively 
large obliquity (Na oz et all 1201 lb . The measured projected 
obliquity is < 16 degrees, implying that either the current true 
obliquity is small, or the star is significantly inclined (i.e., 
/* ~ 0). However, if the star is significantly inclined, then 
the system cannot be synchronized. Perhaps a more likely 
alternative is that, after emplacement by the tertiary and cir- 
cularization of the orbit, the system continued to evolve under 
tidal forces, with KELT-lb migrating inward to its current or- 
bit while damping the obliquity of KELT- 1 and synchronizing 
its spin period. Clearly, detailed simulations are needed to es- 
tablish whether or not this scenario has any basis in physical 
reality. 
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FIG. 18. — Radius versus mass for the known transiting companions 
to main-sequence stars with companion masses in the range 10—100 Mj 
that have measured radii. An estimate of the deuterium burning limit 
I Spiegel et al. 2011) is shown as the vertical dotted line, and the hydrogen 
burning limit is shown as the vertical dashed line. Brown dwarfs are shown 
as triangles, exoplanets as squares, and low-mass stars as asterisks. KELT-lb 
is shown as the large star. Pre dicted radii as a fun ction of mass for isolated 
objects from the isochrones of Baraffe et al. (2003) are shown for an age of 
5 Gyr (dashed), 1 Gyr (dotted), and 0.5 Gyr (long dashed); the true age of 
the KELT- 1 system is almost certainly between 1 — 5 Gyr. Although stellar 
insolation is likely to increase the radii at fixed mass, Bouchy et al. J201 laf) 
predict that the effect is small. KELT-lb therefore has an anomalously large 
radius. 

6.3. Comparison to Theoretical Models of Brown Dwarfs 

Transiting brown dwarfs provide one of the only ways 
to test and calibrate models of BD structure and evolution, 
which are used to interpret observations of the hundreds of 
free floating brown dwarfs for which no direct measurement 
of mass and radius is possible. Given that only 5 transit- 
ing brown dwarfs with radius measurements were previously 
known, KELT-lb potentially provides another important test 
of these models. Figure [18] shows the mass-radius relation 
for the known transiting companions to main-sequence stars 
with companion masses in the range 10-100 My. Being 
close to the minimum in the brown dwarf desert, the mass 
of KELT-lb begins to fill in the dearth of know systems be- 
tween ~ 20-60 Mj up . Furthermore, the formal uncertainty 
in its radius is only ~ 2.5%, thereby allowing for a stringent 
test of models. In contra st, the two transiting BDs with sim- 
ilar masses, CoRoT- 3b dDeleuil et al.l 12008b and KOI-423b 
dBouchy et alJl20lTbh . have much larger radius uncertainties, 
presumably due to the relative faintness of the host stars. 

Evolutionary models for isolated BDs generally predict that 
young (~ 0.5 Gyr) objects in the mass range 10— 100 Mj| lp 
shoul d have radii of ~ Rj up (see the models of iBaraffe et al.l 
2003 in Fig. [T8l . As these objects cool, however, their radii 
decrease, particularly for masses between 50 and 80 Mj up . 
After ~ 1 Gyr, all isolated objects with mass between 20- 
80 Mj up are predicted to have radii < Rj up . The radius we 
measure for KELT-lb is Rp = LIIO^q™ ^Jup, which, at a 
mass of M P = 27.23^ M J up , is ~ 7 a and ~ 10 a larger 
than the radius predicted by Baraf fe et~ai1 (12003b for ages of 



FIG. 19. — Transit depth assuming no limb darkening, i.e., (Rp/R*) 2 , as 
a function of the apparent V magnitude of the host star for a sample of tran- 
siting systems. Brown dwarfs are shown as triangles, exoplanets as squares, 
and low-mass stars as asterisks. KELT-lb is shown as the large star. All 
else being equal, objects in the top left provide the best targets for follow- 
up. KELT-lb has a similar transit depth as the other known transiting brown 
dwarfs, but is significantly brighter. Also labeled are some other benchmark 
systems. KELT-2b (Mp ~ 1.5 Mj llp ) is shown as a large cross (Beatty et al., 
in preparation). 

1 Gyr and 5 Gyr, respectively. KELT-lb is strongly irradi- 
ated, whic h in principle can del ay its cooling and contraction. 
However. lBouchv et al.l d201 lal) predict that the effect of inso- 
lation is small for brown dwarfs in this mass range, although 
their models were for a much more modest insolation cor- 
responding to an equilibrium temperature of 1800 K (versus 
~ 2400K for KELT-lb). Therefore, given the estimated 1.5-2 
Gyr age of the system, KELT-lb is likely to be significantly 
inflated relative to predictions. 

Using the benchmark double tr ansiting BD 2M0535-05, 
Gomez Maque o~Chew et all ((2009) explore models in which 
brown dwarfs have large spots, which reduce the flux from 
their surface, thereby decreasing their effective temperatures 
and increa sing their radii rel ative to those without spots 
(see also Bou chy et al.l d201 lal) ). They find that these can 
lead to significantly inflated radii, but only for large spot 
filling factors of ~ 50%, and for relatively young (~ 0.5 
Gyr) systems. However, a de t ailed spec troscopic analysis of 
that sy stem by lMohanty et alJ (j2010) and Mohant y & Stassuri 
(2012), shows that surface spots cannot be present with 
such a large filling factor, and thus favor global struc- 
tural effects such as strong internal magnetic fields (e.g., 
iMullan & MacDonakl l2010h . Many other mechanisms have 
been invoke d to explain the inflated radi i of some giant exo- 
planets (see iFortnev & Nettelmannll2010l for a review), how- 
ever it is not clear which, if any, of the many mechanisms that 
have proposed may also be applied to inflated brown dwarfs. 

We would be remiss if we did not question whether we were 
erroneously inferring a large radius for the planet. In the past, 
such situations have arisen when there is a discrepancy be- 
tween the constraint on the stellar density from the light curve 
and the constraint on the stellar surface gravity of the star from 
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spectroscopy (e.g JJohns-Krull et al.ll2008tlWinn et"ail 2008). 
In our case, we find no such tension. The parameters of the 
star inferred from the spectroscopic data alone are in nearly 
perfect agreement with the results from the global analysis of 
the light curve, RV data, and spectroscopic constraints. We 
note that the effect of allowing a non-zero eccentricity also 
has a negligible effect on the inferred planetary radius. Fi- 
nally, we reiterate that the faint companion detected in AO 
imaging ( jH3.5K which is unresolved in our follow-up pho- 
tometry, has a negligible effect on our global fit and inferred 
parameters. Therefore, we believe our estimate of Rp is likely 
robust. 

We conclude by noting that there is a need for predictions 
of the radii of brown dwarfs for a range of ages and stellar 
insolations, and it would be worthwhile to explore whether 
or not the inflation mechanisms that have been invented to 
explain anomalously large giant planets might work for much 
more massive and dense objects as well. 

6.4. Prospects for Follow Up 

Figure [19] compares the transit depth and apparent visual 
magnitude of the KELT-1 system (S ~ 0. 6%, V = 10.7) to the 
sample of transiting systems collected in {|6.1.1l with available 
V magnitudes. KELT-1 is not particularly bright compared to 
the bulk of the known transiting exoplanet hosts. However, it 
is significantly brighter than the hosts of all k nown transiting 
brown dwarfs; the next brightest is WASP-30 (I Anderson et al.l 
l20llh . which is - 1.2 magnitudes fainter. On the other hand, 
the depth of the KELT- lb transit is similar to that of the other 
known brown dwarfs. 

The prospects for follow-up of KELT- lb are exciting, not 
only because of the brightness of the host, but also because 
of the extreme nature of the system parameters, in particu- 
lar the relatively short orbital period, relatively large stellar 
radius, and relatively large am ount of stellar irradiation re- 
ceived by the planet. Fo llowing iMazeh & Fa igler (201 (J) and 
iFaigler & Mazehl (1201 lh . we can estimate the amplitudes of 
ellipsoidal variations A e m D , Doppler beaming Abeam (see also 
lLoeb & Ga udi 2003), reflected light eclipses and phase vari- 
ations A re f, and thermal light eclipses and phase variations 

Atherm? 

Abeam = OL b eanA \ — ] ~ 5 .1 a beam X 1 0" 5 (11) 



A re f=a ref ~ 4.6a re/ x 10 4 (12) 

Aeiiip = a em J^- ( ^ j ~4.1a f , //p x 10" 4 (13) 



Afherm ~ Qtherm 



1/2 



3.2a therm x 10" 3 ,(14) 



where the expression for A,h erm assumes observations in the 
Rayleigh-Jeans tail of both objects, and the expression for 
A e iiip assumes an edge-on orbit. The dim ensionless constants 
a are defined in IMazeh & Faiglerl (1201 Ol) . but to make con- 
tact with the secondary eclipse analysis in §!5.4l we note that 
a re f = A g and a,i, erm = [/'(l -As)] 1 ' 4 - All of these constants 
are expected to be of order unity, except for a re f, which may 
be quite lo w for strongly irra diated planets, depending on 
wavelength dBurrows et al.ll2008h . Based on previous results, 
all of these effects with the possible exception of Doppler 
beaming are likely to be detec table with precision photome- 
try (see, e.g. JCowan e t al. 2012). For ellipsoidal variations in 



particular, we expect a e m p ~ 2 and thus a relatively large am- 
plitude of A e iu p ~ 10~ 3 . Furthermore, the detection of all these 
signals is facilitated by the short orbital period for KELT- lb. 

The prospects for transmission spectroscopy are proba- 
bly poorer, given the relatively small planet/star radius ratio 
(~ 0.078) and more importantly the large surface gravity for 
KELT- lb. For the optimistic case of T eq ~ 2400K assuming 
zero albedo and perfect redistribution, the scale height is only 
H r^kT j (pmngp) ~ 16 km, and thus will only lead to changes 
in the transit depth of order ~ 2H/R P ~ 0.04%. 

7. SUMMARY 

We have presented the discovery of KELT- lb, the first 
transiting low-mass companion from the wide-field Kilode- 
gree Extremely Little Telescope-North (KELT-North) transit 
survey. The host star KELT-1 is a mildly evolved, solar- 
metallicity, rapidly-rotating, mid-F star with an age of ~ 
1.5-2 Gyr located at a distance of ~ 260 pc. The transit- 
ing companion is a low-mass brown dwarf or supermassive 
planet with mass ~ 27 Mj up , on a very short period, circular 
orbit of P r~ 1.2 days. 

In many ways, the KELT-1 system is quite unusual and ex- 
treme: KELT- lb receives a large amount of stellar insolation, 
is inflated relative to theoretical predictions, and raises strong 
tides on its host. The obliquity of KELT-1 is consistent with 
zero, and there is evidence that the spin of KELT-1 has been 
synchronized with the orbital period of KELT- 1 . Finally, there 
is a likely M-dwarf stellar companion to the KELT-1 system 
with a projected separation of ~ 150 AU. As the first defini- 
tively inflated transiting brown dwarf, KELT- lb demonstrates 
the need for models of brown dwarfs subject to a range of 
stellar insolations. 

A plausible formation scenario for this system posits that 
KELT- lb formed on a much wider orbit, and was driven 
to a smaller semimajor axis by the tertiary via the Kozai- 
Lidov mechanism. The system then continued to evolve un- 
der strong tidal forces, with KELT- lb migrating inward to its 
current orbit, while damping the obliquity of KELT- 1 and syn- 
chronizing its spin period. 

The future evolution of the KELT-1 system may be spectac- 
ular. As KELT-1 continues to evolve and its radius increases, 
so will the tides raised on it by KELT- lb. Assuming KELT-1 
is and remains tidally locked, as it cools it will develop a deep 
convective envelope, but be forced to rotate at an ever increas- 
ing rate. In ~ 2 Gyr, KELT-1 will have roughly the tempera- 
ture of sun, but with a radius of ~ 2 R and a rotational veloc- 
ity of ~ 100 km s" 1 . At this point, KELT- 1 will l ikely become 
an active RS CVn star (IWalter & Bowverlll98ll) . Eventually, 
as KELT-1 reaches the base of the giant branch, it will swal- 
low KELT- lb wh ole, likely resulting in a bright UV/X-ray and 
optical transient (Metzger et al. 2012b. 
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TABLE 2 
Summary of Observations 



Observatory 



UT Dates 



# Obs Wavelength Section; Table 



Photometry, Primary Transit: 

PvdKO 2011-12-03 — 2011-12-04 

ULMO 2011-12-03 — 2011-12-04 

HAO 2011-12-10 

FLWO 2011-12-16 

ULMO 2011-12-31 —2012-01-01 

FLWO 2012-01-07 

Photometry, Secondary Eclipse: 

ULMO 2011-12-02 
FTN/LCOGT 2011-12-30 
ULMO 2012-01-04 

Spectroscopy, Orbit: 

FLWO/TRES 20 1 1 - 1 1 -09 — 20 1 2-02-05 
Spectroscopy, Rossiter-McLaughlin: 
FLWO/TRES UT 20 1 2-0 1 -07 
Adaptive Optics Imaging: 
Keck/NIRC2 UT 20 1 2-0 1 -07 



151 
110 

266 
362 
162 
459 



115 

72 
121 



23 



PS-Z" 



3900-8900A 



42 3900-8900A 



H, K' 



jt3ATl T[9l 



3.4.2 



3.4.3 



3.4.4 



3.4.2 



3.4.4 



1 10 

in 

Tfl2 

113 
1 14 



3.4.2 



3.4.5 



3.4.2 



EH 10 

jtl3lTf8l 

EH 



NOTE. — ULMO=University of Louisville Moore Observatory; PvdKO=Peter van de Kamp 
Observatory; HAO=Hereford Arizona Observatory; FLWO=Fred Lawrence Whipple Observatory; 
FTN=Faulkes Telescope North; LCOGT=Las Cumbres Observatory Global Telescope Network; 
"Pan-STARRS-Z 



TABLE 3 
KELT- 1 Stellar Properties 



Parameter Description (Units) Value Source Ref. 



Names 




TYC 2785-2130-1 
GSC 02785-02130 
2MASS J000 1269 1+39230 17 






a J2000 




00:01:26.92 


Tycho 


1 


<5j2000 




+39:23:01.7 


Tycho 


1 


Bf 




11.363 ±0.065 


Tycho-2 


2 


V T 




10.701 ±0.057 


Tycho-2 


2 


J 




9.682±0.022 


2MASS 


3 


H 




9.534±0.030 


2MASS 


3 


K 




9.437±0.019 


2MASS 


3 


WISE1 




12.077±0.022 


WISE 


4 


WISE2 




12.732±0.020 


WISE 


4 


WISE3 




14.655±0.030 


WISE 


4 


Ma 


Proper Motion in RA (mas yr~' ) 


-10.1 ±0.7 


NOMAD 


5 


l-l-S 


Proper Motion in Dec. (mas yr~') 


-9.4 ±0.7 


NOMAD 


5 


7abs 


Absolute Systemic RV (km s~') 


-14.2±0.2 


This Paper 






Spectral Type 


F5±l 


This Paper 




d 


Distance (pc) 


262 ±14 


This Paper 






Age (Gyr) 


1.75 ±0.25 


This Paper 




Av 


Visual Extinction 


0.18±0.10 


This Paper 




(U a ,V,W) 


Galactic Space Velocities (km s~') 


(19.9±l.l,-9.6±0.5, -2.6±0.9) 


This Paper 




M c - 


Mass of Convective Zone (Mq ) 


[2.8 ±0.4] x Kr 5 


This Paper 




Rcz 


Base of the Convective Zone (Rq ) 


1.30 ±0.03 


This Paper 




c* 


Total Moment of Inertia (cgs) 


[1.15±0.04] x 10 54 


This Paper 




c cz 


Convection Zone Moment of Inertia (cgs) 


[3.2 ±0.6] x 10 51 


This Paper 





NOTE. — l = H0g et al. 1 1998), 2=Hog et al. (2000), 3= : Skrutskie et al. (2006); Cutri et al. (2003), 4= Wright et al. (2010); 
ICutri & et all 120121) 5= IZacharias et alj)2004l) . 

" We adopt a right-handed coordinate system such that positive U is toward the Galactic center. 
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TABLE 4 

Median values and 68% confidence intervals for the Physical and Orbital 
Parameters of the KELT-1 System 



Variable Description (Units) Value (e =jl 0) Value (e = 0) 



Stellar Parameters: 



M* 


Mass (Mq) 


1.324 ±0.026 


i?* .... 


Radius (R ) 


1 462+° 037 

3 48 +eiJ 
-Q1.S 
f) <q7+0.026 
-Q.Q39 
a 790+0.012 

^■"'-o.ow 




Luminosity (Lq) 


p* 


Density (cgs) 


logg... 


Surface gravity (cgs) 


?eff 


Effective temperature (K) 


6518±50 


[Fe/H]. 


Metallicity 


0.008 ±0.073 


vsin/* . 


Rotational velocity (m s _1 ) 


56000 ±2000 


A 


Spin-orbit alignment (degrees) . . 


2± 16 



Planetary Parameters: 
e Eccentricity . 



ui* Argument of periastron (degrees) 

P Period (days) 

a Semi-major axis (AU) 

M P Mass(Mj) 

R P Radius (Rj) 

pp Density (cgs) 

loggp . . Surface gravity 

T eq ■ ■ ■ ■ Equilibrium temperature (K) .... 

© Safronov number 

(F) Incident flux (10 9 erg s 1 cm 2 ) . 



0.0099+0;™ 
61+71 

1.217514 ± 0.000015 
0.02466 ±0.00016 

27 23 +0 ' 50 

1 110 +0°0 4 3 8 2 
1 - llu -0.022 

24 7+1. f 

A 7<?8+<J^ 17 

-0,023 

2422+g 

g 1 7+0.023 
"•^ -0.028 

7 oi +0.42 

'■ S1 -0.33 



1 077+0.026 
1 452+ ' 033 

\ 45 4 p.w 

0.610+M* 



4 7-JC1+O.UU87 

6517 ±49 
0.009 ±0.073 
56000 ±2000 
1 ± 15 



1.217513 ± 0.000015 
0.02464 ±0.00016 
27.24+°*! 



25.2! I? 
4 744+0:013 
-0.021 

2414+1 
0.919+°? 19 
7.71; 
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TABLE 5 

Median values and 68% confidence intervals for the Lighcurve and Radial 
Velocity Parameters of the KELT-1 System 



Parameter 



Description (Units) 



Value 



RV Parameters: 

Tq Time of inferior conjunction (BJDtdb) 

Tp Time of periastron (BJDtdb) 

K RV semi-amplitude (m s ) 

K R RM amplitude (m s -1 ) 

Mp sin i . . Minimum mass ( Mj ) 

Mp/M, . Mass ratio 

u RM linear limb darkening 

zero point for Orbital RVs (Table[7} (m s _1 ) . 



zero point for RM RVs (Tabled) ( m s_1 ) ■ 



70 — 

71 • • • • 

ecosu* 

esinaj* 

f(m 1 , nil) Mass function ( Mj ) 

Primary Transit Parameters: 

Rp/Rt . . Radius of the planet in stellar radii 

a/R* .... Semi-major axis in stellar radii 

i Inclination (degrees) 

b Impact parameter 

<5 Transit depth 

Tfwhm ■ ■ FWHM duration (days) 

r Ingress/egress duration (days) 

7i4 Total duration (days) 

Pt A priori non-grazing transit probability 

P T q A priori transit probability 

7c o transit time for PvdKO UT 20 1 1 - 1 2-03 (BJDtdb ) 

r cl transit time for ULMO UT 201 1-12-03 (BJDtdb) 

T C2 transit time for HAO UT 20 1 1 - 1 2- 1 (BJDtdb ) ■ 

r C3 transit time forFLWO UT 2011-12-16 (BJDtdb) 

T CA transit time for ULMO UT 201 1-12-31 (BJDtdb) 

r C5 transit time for FLWO UT 2012-01-07 (BJDtdb) 

uisioani ■ ■ linear limb-darkening coefficient 

uzsioani ■ ■ quadratic limb-darkening coefficient 

"isioanr ■ ■ linear limb-darkening coefficient 

u 2Sloam ■ ■ quadratic limb-darkening coefficient 

"isioanz ■ ■ linear limb-darkening coefficient 

u 2Sloanz ■ ■ quadratic limb-darkening coefficient 

Secondary Eclipse Parameters: 

Ts Time of eclipse (BJDtdb) 

bs Impact parameter 

Ts.fwhm ■ FWHM duration (days) 

t$ Ingress/egress duration (days) 

7ji4 .... Total duration (days) 

P$ A priori non-grazing eclipse probability 

Ps.a A priori eclipse probability 



2455914.1628^™23 
2455914.07^'^ 
4239 ±52 
342 ±13 
27 20+ 49 
0.01964 ±0.00028 
U.5S42_„ m9 
-14200 ± 50 (stat.)±200 (sys.) 
-14200+g (stat.)±200 (sys.) 
0.0018+0 ?°p 

-0.00037 



0.0041+0- 
0.01006+r 



0.07801+0 .00060 



3.626+. 080 
87.8+j-3 

I4I+ IUI 

—0 08i 

0.006086+0-±4 

0.10642 ±0.00045 
0.00870;000044 

1151 q+o:ooHAs 
n o««q+o.0074 
2991+0 0088 
2455899.5550 ±0.0010 
2455899.55407 ± 0.00044 
2455905.63860+000084 

245591 1.72553 ±0.00045 
2455927.55574+O-0O040 
2455933 fi43 21+ -000 41 
n 2143+0 00 45 

0.3152S 
0.2856+0-005 

3262+0 0025 
1 64 5+0.8043 

3062+0- 0027 
u - JUDZ -0.0033 



2455913.5560+0-00^ 
0.142+0- 11 
0.1()73+^ 3 

n nns8^+8^8oi7 
u.uu»si_ 00027 

0.1 162+0 0026 

0.2525; 

0.2952 H 



'it 



TABLE 6 
KELT-1 Transit Times 



Epoch 


T C 




O-C 


(0-Q/<tt c 


Obs. 




BJDtdb 


Seconds 


Seconds 






-8 


2455899.55497 


87 


188.12 


2.14 


PvdKO 


-8 


2455899.55407 


37 


110.36 


2.92 


ULMO 


-3 


2455905.63859 


72 


-146.99 


-2.03 


HAO 


2 


2455911.72553 


38 


-195.85 


-5.06 


FLWO 


15 


2455927.55574 


35 


37.16 


1.05 


ULMO 


20 


2455933.64320 


35 


33.84 


0.96 


FLWO 
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TABLE 7 

KELT-1 Orbital Radial Velocities and Bisectors 



BJDxdb 


RV 


CRV 


Bisector 






(m s" 1 ) 


(ms ') 


(ms" 1 ) 


! 1 > 

(ms ') 


94SS874 861 15*5 


-10308 9 


193 3 


103 9 


152 2 


9455876 75115^ 


-18146 


182 9 


121 5 


106 7 


9455884 78^Q08 

ZtJJOOt. / OJ7UO 


1 flQlft 9 


275 9 


1 SQ fi 
1 ov.o 


Qn 1 


9455885 79^564 


- 1 05 1 7 5 


238 3 


56 


140.3 


9455887 750071 
ZtJJoo / . / JUL) / 1 


1 897Q 1 


213 7 


75 6 


130 1 


9455888 70Q8Q1 

itJJ OOO. /U7U71 


-16518 6 

1 \JJ IO.U 


173 4 


-61.1 


117 


9455Q00 658971 


-1 1722 5 


222 


-55 7 


84 4 


9455Q01 77^451 


1 0614 7 




-41 6 


108 3 


9455Q04 764448 


-18416 8 


180 2 


-29.4 


124.5 


9455Q05 5Q8058 
Z4J J VVJJ. J70UJ0 


-131 14 1 






74 8 


9455Q1 1 61 8841 

i'+JJ 7 1 1.U1 Out 1 


-1 1901 7 


240.3 


165.9 


84 3 


2 fill 1 63 

i-t J J7 1 i_ . J 11_J_7 


_i 0950 


325.6 


-33.0 


160.8 


2455930.604440 


-13938.2 


216.2 


67.7 


74.6 


2455931.588925 


-18257.4 


210.1 


131.8 


46.3 


2455932.584867 


-17277.9 


190.9 


-158.3 


75.2 


2455934.572100 


-10043.3 


248.2 


-53.4 


82.5 


2455935.671466 


-11237.1 


334.9 


-293.1 


130.6 


2455936.607043 


-16406.3 


307.3 


-34.7 


209.3 


2455937.570378 


-18299.2 


156.2 


4.3 


68.5 


2455940.565265 


-10308.6 


224.9 


-210.3 


102.9 


2455957.622238 


-10008.5 


267.8 


-123.9 


97.2 


2455960.657124 


-17869.8 


297.0 


10.2 


171.4 


2455962.599920 


-9972.1 


235.1 


9.8 


121.1 



NOTE. — Radial velocities and bisector spans with uncertainties for 
spectra taken with TRES at FLWO on nights outsi de of transit. Times ar e 
for mid-exposure and are in the BJDtdb standard 1 Eastman et al. 2010b). 
RVs are approximately on an absolute scale, and the uncertainties have 
been scaled based on a global fit to the data. See fl3.3| and fl5. 2 1 for 
details. 
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TABLE 8 

KELT- 1 RM Radial Velocities and Bisectors 



BJDxdb 


RV 


&RV 


Bisector 






(m s" 1 ) 


(ms l ) 


(m s" 1 ) 


(ms [ ) 


94.SSQH Sfi14.fiS 


-12405 3 


181.6 


53 5 


109 8 


Z,tJJ 7JJ .JUJ Itu 


-12596.4 


217 


-178 


141 5 


94.SSQH SfiSS^S 


-12713 2 


901 o 

ZWj.U 


10 2 


145 3 


94.SSQH S719R9 


-12618.2 


191 4 


218 3 


130 1 


94.SSQH S7fiQ71 


-12958 8 


160 9 


01 7 


122 5 


94.SSQH ssnfisd 


-13001 


215 5 


11 1 


151 7 


94.SSQH SRS917 


-12937 2 


169.0 


-78 


124 5 


94.SSQH SSRQ17 


-12761 6 


153 6 


-21 1 


14.Q i 


94.SSQT3 SQ9fiOQ 


-12974 8 


181 8 


-59.1 


108 3 


94.SSQH SQfiQRl 


1 9Q4.S s 


1 SS 8 

1 JJ.O 


so 


1 14. s 


94.SSQH ftOOf^ftd 

ZHJJ J J J .uuuuut 


-12643 7 


112 7 


-74 2 


81 3 


94.SSQH fi04.1SS 

ZtJJ 7JJ .UUt'JJJ 


-131 10 1 


191 6 


-34 8 


113 6 


94.SSQH fiOS7fiS 
Z4JJ7JJ.OUO / UJ 


1 HOfi s 

- 1 jjUD.O 


Q4. 


9fifi Q 


129 1 


94.SSQH fi1 94.S7 

ZtJJ 7JJ.U1 Zt J / 


-1291 1 2 


136 3 


-4.4 


122 2 


94.SS0^^ filfil 17 

ZtJJ 7JJ.U1U1J / 


-13378 9 


112 5 


90.8 


63 


94.SSQ11 fi91 104. 

ZtJJ7jJ.UZl 171 


-13577 2 


153 6 


148 1 


114 1 


ZtJJ 7 J J . UitO / J 


-13770.6 


129.3 


53.1 


99.4 


94.SSQ11 fi9SSSS 
Zt J J V J J .DZOJJ J 


1 1R0Q 


140 7 


1Q s 

17.J 


1 17 1 


94.SS0^^ fv*3fR4 

ZtJJ 7JJ . U J JV7 Jt 


- 1 J7JJ.U 


139 4 


22 5 


115 4 


94.SSQH filfi714. 

ZtJJ 7JJ.UJU/ It 


-13897 2 


115 7 


177 2 


130.4 


94.SSQ11 fi4.04.Ofi 

ZtJJ 7 J J .Utl/tlJU 


-14256 


132 6 


45.2 


124 7 


94.SS0^^ fidSfHS 

ZtJJ 7JJ .UtJV/JJ 


-14345 2 


143 2 


0.4 


82 1 


94.SSQ11 fi4.S7 1 S 


14.179 Q 
-It J 1 L.y 


145 8 


90 Q 

■ZI/.7 


88 2 


94.SSQH fiS91Qfi 

ZtJJ J J J .UJZJ7U 


-14371 8 


173 9 


-188 


86.4 


94.SSQH fiSfiSfiQ 

ZtJJ7JJ.UJUOU7 


-14573 1 


148 5 


10 8 


87 1 


94.SSQH fifiOSSS 

ZtJJ 7JJ . UU\/J J J 


-14845 4 


137 6 


-234 2 


104.5 


94.SSQH fifi4.91S 

ZtJJ 7JJ .UUtZJJ 


-14862 9 


126 8 


-117 8 


71 


94.SSQ11 ftf\RHRl 

ZtJJVJJ . UUOOO / 


14.Q7fi l 

-it? / U. J 


136 4 


1 Si fi 

- 1 oj.U 


88 4 


94.SSQT3 fi79Sfi7 

ZtJJ 7JJ.U / Z iyJ I 


-15112 1 


173 7 


37 9 


136.2 


94.SSQH fi7fi94.S 

ZtJJ 7JJ.U / uzto 


-1 S071 S 
- 1 j i j / j . j 


104.0 


-21.6 


112 4 


94.SSQH fiS07fi1 

ZtJJ 7JJ.UOU/U1 


-15482 9 


140 4 


2118 


76.4 


2455933.684453 


-15385.5 


181.2 


-82.2 


198.6 


2455933.688133 


-15349.1 


141.5 


52.1 


108.0 


2455933.692820 


-15614.2 


156.9 


164.2 


135.4 


2455933.696501 


-15427.2 


206.4 


-237.5 


142.5 


2455933.700192 


-15335.8 


313.0 


-699.1 


412.4 


2455933.704648 


-15667.5 


145.1 


157.8 


119.1 


2455933.708340 


-15366.4 


218.3 


-54.8 


133.8 


2455933.712020 


-15815.1 


219.6 


181.9 


114.5 


2455933.716464 


-15824.6 


220.1 


-206.9 


147.1 


2455933.720145 


-16127.7 


180.2 


-245.7 


210.6 


2455933.723836 


-15832.2 


182.2 


-33.8 


160.4 



NOTE. — Radial velocities and bisector spans with uncertainties for 
spectra taken with TRES at FLWO on the night of the primary transit 
on UT 2 012-01-07. Times are f or mid-exposure and are in the BJDtdb 
standard i Eastman et al. 2010b). RVs are approximately on an absolute 
scale, and the uncertainties have been scaled based on a global fit to the 
data. See g|l and ^2|for details. 



KELT- lb 



TABLE 9 

Relative Photometry from PvdKO on 
UT 201 1-12-03 (0 



BJDtdb 


Relative Flux 


Uncertainty 


2455899.454080 


1.00112 


0.00177 


2455899.455342 


0.99986 


0.00177 


2455899.456603 


0.99606 


0.00206 


2455899.457865 


0.99556 


0.00218 


2455899.459138 


1.00147 


0.00192 


2455899.460399 


0.99946 


0.00163 


2455899.461661 


0.99686 


0.00162 


2455899.462934 


1.00004 


0.00190 


2455899.464207 


1.00294 


0.00221 


2455899.465468 


0.99780 


0.00175 



NOTE. — This photometry has been corrected for a 
linear trend with airmass, and normalized by the fitted 
out-of- transit flux. See j j5.2| This table is published in its 
entirety in the electronic edition. A portion is shown here 
for guidance regarding its form and content. 



TABLE 10 

Relative Photometry from ULMO on UT 
2011-12-03(r) 



BJDtdb 


Relative Flux 


Uncertainty 


2455899.475596 


1.00155 


0.00091 


2455899.476984 


1.00058 


0.00091 


2455899.478361 


0.99852 


0.00089 


2455899.479738 


1.00015 


0.00091 


2455899.481115 


1.00175 


0.00089 


2455899.482493 


1.00039 


0.00089 


2455899.483870 


1.00020 


0.00089 


2455899.485259 


1.00144 


0.00089 


2455899.486636 


0.99967 


0.00089 


2455899.488014 


1.00099 


0.00091 



NOTE. — This photometry has been corrected for a 
linear trend with airmass, and normalized by the fitted 
out-of-transit flux. See i |5.2| This table is published in its 
entirety in the electronic edition. A portion is shown here 
for guidance regarding its form and content. 



TABLE 11 

Relative Photometry from HAO on UT 
2011-12-10(0 



BJDtdb 


Relative Flux 


Uncertainty 


2455905.540572 


0.99833 


0.00250 


2455905.541132 


0.99667 


0.00250 


2455905.541712 


1.00321 


0.00252 


2455905.542272 


1.00209 


0.00243 


2455905.542842 


1.00116 


0.00243 


2455905.543412 


0.99931 


0.00243 


2455905.543972 


1.00078 


0.00243 


2455905.544542 


0.99829 


0.00242 


2455905.545112 


0.99994 


0.00243 


2455905.545692 


0.99800 


0.00250 



NOTE. — This photometry has been corrected for a 
linear trend with airmass, and normalized by the fitted 
out-of-transit flux. See fl5.2| This table is published in its 
entirety in the electronic edition. A portion is shown here 
for guidance regarding its form and content. 



Siverd et al. 



TABLE 12 

Relative Photometry from FLWO on UT 
201 1-12-16 (z) 



BJDtdb 


Relative Flux 


Uncertainty 


2455911.612185 


1.00109 


0.00154 


2455911.612983 


1.00230 


0.00154 


2455911.613516 


1.00009 


0.00154 


2455911.614002 


1.00177 


0.00154 


2455911.614488 


0.99899 


0.00154 


2455911.615078 


1.00118 


0.00154 


2455911.615564 


1.00159 


0.00154 


2455911.616051 


0.99953 


0.00154 


2455911.616560 


1.00155 


0.00154 


2455911.617081 


1.00066 


0.00154 



NOTE. — This photometry has been corrected for a 
linear trend with airmass, and normalized by the fitted 
out-of- transit flux. See j j5.2| This table is published in its 
entirety in the electronic edition. A portion is shown here 
for guidance regarding its form and content. 



TABLE 13 

Relative Photometry from ULMO on UT 
2011-12-31 (r) 



BJDtdb 


Relative Flux 


Uncertainty 


2455927.468704 


0.99844 


0.00145 


2455927.470335 


0.99820 


0.00125 


2455927.471713 


0.99925 


0.00116 


2455927.473090 


1.00005 


0.00105 


2455927.474467 


0.99994 


0.00099 


2455927.475844 


1.00105 


0.00092 


2455927.477233 


0.99940 


0.00090 


2455927.478610 


1.00084 


0.00086 


2455927.479987 


0.99960 


0.00086 


2455927.481365 


0.99849 


0.00084 



NOTE. — This photometry has been corrected for a 
linear trend with airmass, and normalized by the fitted 
out-of-transit flux. See i |5.2| This table is published in its 
entirety in the electronic edition. A portion is shown here 
for guidance regarding its form and content. 



TABLE 14 

Relative Photometry from FLWO on UT 
2012-01-07 (0 



BJDtdb 


Relative Flux 


Uncertainty 


2455933.568776 


1.00164 


0.00184 


2455933.569101 


0.99893 


0.00184 


2455933.569459 


1.00032 


0.00184 


2455933.569830 


1.00046 


0.00184 


2455933.570200 


0.99947 


0.00184 


2455933.570536 


1.00032 


0.00184 


2455933.570871 


1.00131 


0.00184 


2455933.571195 


1.00074 


0.00184 


2455933.571531 


0.99764 


0.00184 


2455933.571867 


1.00003 


0.00184 



NOTE. — This photometry has been corrected for a 
linear trend with airmass, and normalized by the fitted 
out-of-transit flux. See fl5.2| This table is published in its 
entirety in the electronic edition. A portion is shown here 
for guidance regarding its form and content. 



KELT- lb 



TABLE 15 
Relative Photometry from ULMO on 
UT 2011-12-02 (i) 



BJDtdb 


Relative Flux 


Uncertainty 


2455897.65499 


1.00097 


0.00145 


2455897.65721 


0.99944 


0.00134 


2455897.65859 


0.99993 


0.00140 


2455897.65997 


0.99930 


0.00141 


2455897.66134 


1.00078 


0.00140 


2455897.66273 


1.00023 


0.00138 


2455897.66411 


1.00137 


0.00134 


2455897.66549 


0.99971 


0.00138 


2455897.66686 


0.99978 


0.00140 


2455897.66824 


0.99970 


0.00140 



NOTE. — This photometry has been normalized and 
corrected for a linear trend with time. See fl5.4J This 
table is published in its entirety in the electronic edition. 
A portion is shown here for guidance regarding its form 
and content. 



TABLE 16 

Relative Photometry from FTN/LCOGT 
on UT 201 1-12-30 (PS-Z) 



BJDtdb 


Relative Flux 


Uncertainty 


2455925.730250 


0.99855 


0.00078 


2455925.731230 


1.00045 


0.00078 


2455925.732210 


0.99953 


0.00078 


2455925.733190 


1.00032 


0.00078 


2455925.734190 


1.00069 


0.00078 


2455925.735170 


1.00000 


0.00078 


2455925.739770 


1.00089 


0.00078 


2455925.740700 


0.99925 


0.00078 


2455925.741670 


0.99874 


0.00078 


2455925.742660 


0.99933 


0.00078 



NOTE. — This photometiy has been normalized and 
corrected for a linear trend with time. See j j5.4J This 
table is published in its entirety in the electronic edition. 
A portion is shown here for guidance regarding its form 
and content. 



TABLE 17 
Relative Photometry from ULMO on 
UT 2012-01-04 (/) 



BJDtdb 


Relative Flux 


Uncertainty 


2455930.50268 


1.00020 


0.001 1 1 


2455930.50487 


1.00126 


0.00110 


2455930.50624 


0.99911 


0.00110 


2455930.50762 


0.99970 


0.00111 


2455930.50901 


1.00139 


0.00116 


2455930.51039 


1.00146 


0.00111 


2455930.51177 


0.99958 


0.00111 


2455930.51314 


1.00016 


0.00111 


2455930.51452 


1.00098 


0.00111 


2455930.51590 


1.00085 


0.00111 



NOTE. — This photometry has been normalized and 
corrected for a linear trend with time. See jj5 .4 \ This 
table is published in its entirety in the electronic edition. 
A portion is shown here for guidance regarding its form 
and content. 



